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ABSTRACT 


This thesis is concerned with the development of a model of 
an eight tray distillation column and subsequent use of tnis model for 
simulation and experimental testing of a multivariable control system. 

The first portion of this work is a literature survey of 
distillation modelling techniques and of distillation column control 
strategies. A description and a discussion of many of the more common 
assumptions used in modelling are presented. The models appearing in 
the technical literature are considered with special attention being 
given to the assumptions that have been utilized and application of 
the models to control studies. Many of the control strategies which 
appear in the technical literature are presented, compared and evalu- 
ated on the basis of the necessity of a model and also if any experi- 
mental verification exists. 

The twentieth-order linear state space model of the pilot 
scale distillation column is derived by using the transient mass and 
energy balances about each stage in linear perturbation form. The 
changes in internal liquid flows are included as a first order function 
of time. The states in this model are the stage enthalpies and the 
liquid flows from each stage. Experimental evaluation indicates good 
agreement between the behavior of the model and the pilot scale equip- 
ment. Two model reduction techniques are presented and applied to the 
model to obtain a low order model for control studies. A fourth-order 
and a second-order model are obtained and the latter is selected for 
the control application. 

Multiloop control studies indicate that the selected product 


control loops exhibit interactive effects, and so a multivariable control 
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approach is applicable. The second-order model, with the states being 
the product enthalpies is used to derive a multivariable control law. 

It is found that measurement noise and time delays in the bottoms com- 
position measurement from the gas chromatograph system cause instability 
in tne bottom enthalpy control loop. Simulation studies on an adjusted 
feedback and a feedback-feedforward control system indicate moderately 
good controlled behavior, with the bottoms composition exhibiting sus- 
tained oscillations. Subsequent experimental runs verify the appli- 
cability of the control schemes and the results predicted by the 


Simulation studies. 
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CHAPTER I 
INTRODUCTION 


The approach to control systems most extensively used in the 
process industries nas been through the use of pneumatic or electronic 
controllers in a multiloop configuration. This approach nas been used 
despite the problems of control variable selection and control loop 
interaction because of the inherent complexity of many processes. This 
complexity has been coupled with the liberal design of processing equip- 
ment to remove the need for exacting control systems. However, the 
recent increase in the number of industrial real time computers and 
advances in process modelling have increased the feasibility of multi- 
variable control system applications. 

The usual objective of a control system applied to a distilla- 
tion column is to maintain product compositions within specific limits 
as dictated by downstream processes or required product specifications, 
despite input disturbances. In the case of a column with no sidestream 
drawoff a typical multiloop product control system could include: 

i) Bottom product composition control by manipulation of 

heat input to the reboiler. 

ii) Top product composition control by manipulation of reflux 

flowrate. 
However, the simultaneous application of both control loops can lead to a 
deterioration in control of both compositions due to the interaction 
between the two loops. This interaction arises from the effects of the 
manipulated variable of one loop on the controlled variable of the other 
loop. The presence of this interaction indicates that the control of a 


distillation column is an ideal situation for the application of a 
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multivariable control scheme. 

However, for the design of a multivariable control system, it 
is necessary to be able to predict the interactive effects. This, in 
general, requires a process model which includes the effects of all 
disturbances as well as these interactive effects. Two general approaches 
to the dynamic modelling of distillation columns which have been found 
successful are: 

i) The application of the governing differential equations 

to yield a model comprised of a system of coupled first 
order differential equations. 

ii) The application of experimental transient data to yield 
an approximate transfer function model. 

For the pilot-scale distillation column located in the Depart- 
ment of Chemical and Petroleum Engineering at the University of Alberta, 
the latter approach has been successfully applied to modelling and to 
the subsequent design of control systems [56,57]. The objective of 
this study will be to use the former approach to obtain a process model, 
and to apply that model to the design of a multivariable contro] system. 
Tne distillation column referred to above will be used in this study to 


experimentally evaluate the control system. 
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CHAPTER II 
LITERATURE SURVEY 


2.1 Introduction 

Tne control of the compositions of products from a distilla- 
tion column has been shown to be quite difficult. Rosenbrock [42] has 
presented a general discussion of problems encountered, and Davison [13] 
has extensively studied one of these problems, that of control loop 
interaction. The development of control strategies for distillation 
Systems has received extensive discussion in the technical literature. 
The more recent design of control strategies has hinged on the develop- 
ment and application of a process model, and so this has also received 
attention in the technical literature. This chapter will attempt to 
outline several approaches to modelling and subsequent control which 


nave been applied. 


2.2 Survey of Distillation Modelling Techniques 

In the technical literature, two general] approaches to distil- 
lation modelling can be distinguished. The first of these methods 
involves the application of the dynamic mass, component and energy 
balances which apply to sections of the system. These equations are 
then simplified by the application of certain assumptions and by the 
knowledge of specific physical property data for the system involved. 
The second approach to modelling involves the use of experimental data 
from a specific system to obtain an approximate transfer function 
relationship between the relevant variables. In this section, several 
models developed from the two different approaches will be compared and 


discussed with specific attention being paid to application and model 
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testing. 

In an early series of articles, Rosenbrock [39-41] has presented 
the general form of the transient balance relationships. He also dis- 
cusses, at some length, several simplifying assumptions, metnods of 
solving the resulting equations, and the development of computer pro- 
grams for this purpose. Subsequent to this series of articles, review 
articles by Archer and Rothfus [3] in 1961 and by Williams [49] in 1963 
have discussed the progress that has been made and have given an evalua- 
tion of the specific modelling approaches proposed. 

In recent years, several theoretically based models, established 
using various different assumptions have been presented in the literature. 
A brief description of some typical assumptions and some remarks on their 
Significance are presented in Table 2.2-]. A number of models which have 
been proposed are classified in Table 2.2-2, with some brief remarks on 
the type of model, its application and specific assumptions made in the 
derivation. It should be noted that in the case of results of specific 
groups of workers, for example the University of Delaware and the Univer- 
sity of Florida groups, several papers have been presented together. 

Svrcek [46] has presented a complicated model solution, which 
shows very good agreement with experimental data. In many cases [4,5,45, 
26 ,20,50,51,52] the binary system used was chosen to enable several 
Simplifying assumptions to be made. The most common choice is two com- 
ponents with similar heats of vapourization, allowing the assumption of 
equimolal overflow to be made. Davison [12] has attempted to define the 
proolem in terms of pressure variations but many simplifying assumptions 
and approximations were necessary in order to obtain a feasible model. 


Several other authors have proposed similar models [1,22], but those 
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presented in Table 2.2-2 represent the best documented and tested models. 
The other approach to distillation column modelling has been 
through the use of experimental data, and often the use of pulse test 
data. Janis [24] and Wood and Pacey [56] have used pulse data, analyzed 
by the use of Fourier transforms to arrive at a Bode plot. An arbitrary 
form of transfer function is then employed to represent the data. Wood 
and Pacey [56] indicate good agreement between experimental transients 
and transfer functions of the first order plus time delay form. Chanh 
[10] has presented and compared several techniques for the time and fre- 
quency domain fitting of experimental data. He and Berry [6] have applied 
time domain methods and indicate good agreement with experimental data 


for simple transfer function forms. 


2.3 Survey of Model Reduction Techniques 

The problem of reducing a state space model involves the removal 
of specific states because they are unmeasurable or don't contribute 
greatly to the overall transient solution. Several methods have been 
proposed in the literature to accomplish this end. Marshall [31] and 
Davison and Chadha [14] have each presented methods based on a modal 
analysis of the system. The difference between the two methods is in the 
approximation used to account for the dynamics of the states to be 
removed. Wilson et. al. [53] have discussed these methods at some length, 
and extended each case to the discrete form of the model. Anderson [2] 
has presented a different method which is actually a least squares 
regression fit of the desired states. These methods will be discussed 
and compared further in the chapter of this work dealing specifically 


with this topic. 
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2.4 Survey of Distillation Column Control Methods 

Historically distillation column control problems have been 
simplified by the use of an over-designing policy for process equipment 
so tnat utilization of simple control systems gave adequate results. 
However, with the increase in knowledge of distillation system behavior, 
and the development of on-line analysis techniques and control systems, 
more stringent design and control methods are being applied. 

Because of the problems which occur in trying to measure 
compositions on a continuing basis, the use of temperature as a composi- 
tion indicator has been widely used. Its application is based on the 
assumption that a mixture composition can change only slightly without 
affecting the mixture temperature, even in a multicomponent mixture, 
where this is not necessarily the case. Williams [50] has outlined the 
estimation of compositions from temperature measurements in a binary 
system, and presents some problems encountered in its application to 
control. The use of product temperature as the controlled variable has 
been shown to yield poor control [47]. In order to minimize product 
composition changes, an intermediate tray temperature has been suggested 
as a controlled variable. An early analysis by Boyd [8,9] indicates 
several possible control schemes utilizing intermediate tray temperatures. 
Similar analyses have been presented by Tivy [47] and by Pyle [34]. Al] 
authors indicate that the temperature sensor used in a control loop 
Should be placed between the feed tray and the end of the column where 
the composition is being controlled. Furthermore the sensor should be 
placed on the tray where a steady state analysis of the column temperature 
gradient, indicates the largest slope. However, the sensor should also 


be as close to the end of the column as possible in order to minimize the 
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time lag between the measured temperature change and the expected product 
composition change. Wood [55] follows a similar approach to sensor 
location but extends the problem to include pressure control by control- 
ling the temperature difference between the top tray and an intermediate 
tray. Chanh [10] has studied extensively the problem of sensor location 
with special consideration given to sensitivity and steady state error. 
He presents simulation results from two distillation systems and subse- 
quent experimental results from one of these systems. In general, the - 
approach of using an intermediate tray temperature has the advantage 
that a control systen can be designed without a process model and using 
only steady state temperature data. However, many problems can arise as 
indicated by several of the above authors. 

As indicated in Section 2.1, the major problem encountered in 
a two point control problem, that is, simultaneous control of both end 
product compositions, is the interaction between control loops. Several 
authors, Luyben [28] and Zalkind [59,60] have proposed methods to com- 
pensate for this interaction. Niederlinski [33] has presented an analysis 
and a criterion to indicate whether this decoupling is advantageous. The 
decoupling method entails the use of a process model to design "inter- 
action compensators" for the control loops. Berry [6] has applied this 
approach to the design of a noninteracting control system for the distil- 
lation column used in this study, using a transfer function model derived 
from pulse test data. He reports a vast improvement from the simple 
multiloop interacting case to the noninteracting case. Berry [6] has 
also applied a method suggested by Rijnsdorp and van Kampen [36-38] which 
nas a decoupling effect. The method is based on controlling the ratio 


of the reflux flow rate to overhead vapour flow rate. The actual ratio 
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used is varied according to measured composition changes. This ratio 
control tends to remove the effect of boil-up rate on the overhead 
control loop thus removing part of the interaction. Berry reports that 
this control scheme also gave improved control performance, being 
approximately the same as the noninteracting case. This method also 
has an advantage over the previous method in that no process model is 
required for compensator design. 

The control schemes discussed thus far have been, in general, 
some form of feedback control, that is correcting for changes in the 
controlled variable and not accounting for any disturbances which may 
enter the system. This latter action is designated as feedforward 
compensation, and its inclusion usually improves control performance 
Since no change in controlled variable need occur before some corrective 
action is taken. However, it has the added requirement that some know- 
ledge of the effect of each disturbance on the behavior of the system is 
needed. One method of accomplishing feedforward control which has been 
Suggested by Luyben [27] is by manipulation of the feed inlet location. 
A subsequent experimental study has been reported by Speicher and 
Luyben [44]. These workers used feed plate manipulation to account for 
feed composition changes. The steady state effects of feed plate loca- 
tion and of feed composition were investigated and this information used 
in the design of a control strategy. They report moderately good control 
results, however, it was found necessary to introduce a time lag into the 
control action. They also point out that the option of using several 
different trays as the feed tray is required in order to prevent over- 
compensation. 


Distefano [16] describes the modelling and subsequent feedforward 
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control of a distillation column, with overhead composition being the 
controlled variable. His model [15,16] was modified slightly to yield 

a form which could be used for control studies. The reflux flow rate 

was used as a control variable to compensate for feed composition 
changes. The required reflux flow rate as a function of time was cal- 
culated beforehand for a given feed composition change, and then manually 
set during a corresponding experimental test. He reports good control of 
the system but some discrepancy between the simulated and experimental 
controlled transients, which he attridutes to reboiler modelling errors. 

Janis [24] has conducted a brief study of feedforward control- 
lers derived from different models of a distillation column. He used a 
model derived from experimental pulse data and a model from linearized 
perturbation equations. He compares the models and controllers but 
presents no experimental evaluation of the control] systems. The use of 
a feedback-feedforward control law derived from a transfer function model 
of a column has been described by Wood and Pacey [56]. The model was 
derived from pulse data. The data were fit in the frequency domain to 
a first order plus time delay function. The controllers derived from 
these transfer functions were used to control overhead product composi- 
tion when the column was subject to feed flow disturbances. The reported 
experimental results indicate good control performance, when feedforward 
action is included. 

A great deal of effort has recently been spent on the control 
and optimization of systems of the linear time invariant state space 
form. Newell [32] has outlined to some extent the derivation of such 
optimal control laws and their application to the process industry. 


Brasilow and Handley [9] have presented the optimal feedback case for a 
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quadratic performance index as the solution to the matrix Ricatti equa- 
tion. The model used was obtained by linearizing the balance equations 
about the steady state operating point. The experimental tests presented 
used temperature as a composition indicator, and the controlled transients 
demonstrate good control and only small offsets for step load changes. 

Hu [22] has presented an extensive review of control strategies for 
distillation columns. He has applied dynamic programming to the optimiza- 
tion of a linearized balance model and presents simulation results to 


demonstrate and compare the control laws. 
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CHAPTER III 


MODEL DEVELOPMENT 


3.1 Introduction 

In order to obtain a model of any process which can be mathe- 
matically solved, it is usually necessary to make simplifying assump- 
tions, which hopefully can be justified by the physical system, and 
also will not restrict the solution. The assumptions which will be 
applied in the development of the distillation column model used in 
this study are: 

1) The vapour dynamics are negligible compared to the liquid 

dynamics. 

2) For small changes in steam rate to the reboiler, the 

changes in vapour flows are negligible. 

3) The liquid and vapour phases in each stage are each well 

mixed, and so are homogeneous in composition and 
temperature. 
) The time rate of change of vapour mass holdup is small. 
5) The heat losses from each stage are constant with time. 
) The pressure throughout the column is constant. 

7) The fluctuations in time dependent variables are small 

so that linearized perturbations can be used. 

The above assumptions are based either on the physical system 
or on tne observations of other workers. The sixth assumption is 
justified by the presence of a pressure controller on the condenser 
using the cooling water flow as the control variable. The last assump- 


tion allows the time varying properties to be expressed in the form: 
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G = G+G' 


where G = any general time varying property 


(ep) 
I 


the "steady state" value of G 

G' = the deviation in G from the "steady state" value 
Since the perturbations are assumed to be small, the terms involving 
products of perturbation terms can be neglected. 

A schematic diagram of the distillation column being considered 
in this study is shown in Figure 3.1-1. The column has eight bubble-cap 
trays, a basket-type steam reboiler, and a total condenser. The feed is 
introduced onto the fourth tray under flow and temperature control. The 
product streams are drawn off under level control on the appropriate 
vessel. The dynamic equations describing each stage are outlined in the 


following sections, and the nomenclature used is presented in Appendix A. 


3.2 Tray Model 

The trays of a distillation column represent a group of stages 
which have the same set of governing equations and so a general set of 
equations can be derived. Figure 3.2-1 shows a schematic representation 
of a tray, designated as stage I. The unsteady state total mass and 


energy balances on this general stage yield: 
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These expressions can be transformed by applying the steady state rela- 


tions, and the perturbation relations of the form: 
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PIGURE S.1-1: 


SCHEMATIC REPRESENTATION OF DISTILLATION COLUMN 
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(3.2-2) 


These equations apply to each tray in the distillation column except the 


feed tray, where terms 


the feed stream. 
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The equations which apply to the feed tray (stage 5) 
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FIGURE 3.2-1: SCHEMATIC REPRESENTATION OF A TRAY 


2] 


3.3 Condenser Model 

The condenser on the equipment under study is the total con- 
denser shown schematically in Figure 3.3-1. The heat removed from the 
condenser by the cooling water is designated as Q and is a load vari- 


able in the model. The total mass and energy balances yield: 
d 
dt “Lio * Nero! 


VQH 


d ; ; : 
dt (yoPio * Weqotio) = Vol - Ly9y9 - Dh 


10 ~ % ~ %o 
Substitution of the perturbation expressions and the steady state mass 


and energy balance relationships yields the following equations: 


Monat Wen ie loo lng oD (3.3-1) 
Werotig * HioWe10 * SLio10 * yo"LI0 
= Vgty + HoVg - Ayollyg + D) - hyg(lyy + 0) - Q (3.3-2) 


3.4 Reboiler Model 
The reboiler used in this study is a basket type reboiler, shown 
schematically in Figure 3.4-1. The mass and energy balances about this 


Stage yield the following expressions: 
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Application of the linearized perturbation and steady state relations 


yields: 
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FIGURE 3.3-1: SCHEMATIC REPRESENTATION OF THE TOTAL CONDENSER 


/ E +i = 7 a’ 
\ y, , 
Se 
| / \ 
f \. 
bi es 
a rr 
th . — 
—___——_}+-——> | 
ya \A/ j 
fs ‘, v¥ 
C3) On nae 
——— v4 Ps 
{ : 5 VV j 
: \ Of Jf 
| ‘ 
~~ i> 
i ; 
\ t 
-" ; 
7. Wy 
2 | cael ———¥ 
—— a - 
7 
7 7 : 


ss : : A ,. » 7 
weaeusding JATOT ANT FO WOUTATHARIAGIN DUTANIMDe .~80e- anUONy 


—— 
<> + 
7 


Aa 


EIGURE S4—=10 


ae 


SCHEMATIC REPRESENTATION OF THE REBOILER 
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WM poe be ee (3.4-1) 
Weqmy + PW, + Wop Hy + HyWey 
= Lofty + holy ~ VyHy - HyV, - Cyhy = hyL; + Q, (3-4-2) 


3.5 Matrix-Vector Representation 

The linearized perturbation balance expressions presented in 
the previous sections represent two sets of ten simultaneous differen- 
tial equations, one set derived from the total mass balance, tie other 
from the energy balances. These expressions can be written as two 
matrix-vector equations as outlined in detail in Appendix B. These 


equations are: 


D (3.5-1) 
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where Wo» W >» H, h, te V are diagonal matrices of steady state values; 


WS 


; B, Cs Co are appropriate coefficient matrices, and the vectors con- 
Sist of perturbation terms. An analysis of these equations indicate 
that there are 20 scalar relationships in 60 time dependent unknowns. 
If it is assumed that: 

1) the vapour dynamics are negligible compared to the liquid 


dynamics, 
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2) for small changes in steam rate to the reboiler, the 
changes in vapour flow are negligible, 
then the following simplifications can be applied. 
i) Ne = Q 
dae 0 


Substitution of these approximations into Equations (3.5-1) and (3.5-2) 


and the subsequent elimination of the vector aH yields: 


WH +h = ALD+CAH-HAL+BVH+(C,-hCD  — (3.8-3) 


By applying these assumptions and substitutions the problem has been 
reduced to one of ten equations in 30 unknowns. Therefore, in order to 
solye the system of equations, 20 auxiliary equations are required. 
Since H and h are the vectors of enthalpies in each stage, it can be 
assumed that they are related by some form of pseudo-equilibrium 
relationship. The vector L is composed of the liquid flow perturba- 
tions, and so should be related in some way to the input liquid flow 
disturbances. These two simplifying assumptions are pursued in the 
following sections to yield the relationships required to allow solution 


of the equations. 


3.6 Enthalpy Linearization 

Since in a binary system, with pressure constant, or specified, 
there is only one degree of freedom, it seems reasonable to assume that 
the vapour and liquid enthalpies should be related. If the assumption 
tnat the "tray efficiency" is constant is made, then it is possible to 


express a linearized enthalpy relationship as: 
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H = by nes b. 


From this, the perturbation form can be expressed as: 


H = by h 


This can be shown as follows: 


In general H = f,(y,T) 
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Furthermore, the assumption of a constant efficiency allows the use of 


‘pseudo-equilibrium' relations of the form: 


So that the enthalpies can be expressed as: 
H = fely) 
= F(x) 


For small perturbations in the variables, the term by can be expressed 


as: 


If a pseudo-equiliorium relationship of the form: 
y = F5(x) 


1s used, then by can be expressed as: 
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This term, bi» the enthalpy linearization factor, can be evaluated for 
eacn stage, knowing the steady state values and the form of the enthalpy 
and ‘equilibrium' functions. Sample calculations of this factor are 
presented in Appendix C. When all stages are being considered, the 


enthalpy linearization can be written as a vector-matrix form: 
H = 8B, h (3.6-1 
H = Bh ) 


where 1 is a diagonal matrix of the enthalpy linearization factors. 


— 
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3.7 Liquid Flow Dynamics | 

Several authors [40-42, 3-5, 25,26] have stated that the 
dynamics of the liquid flow rates within a distillation column can be 
Significant. In general, the liquid flowrate from a tray is dependent 
on tne tray design, the liquid and vapour flows to that tray, and the 
fluid properties. As several workers have outlined [4,5,25,26], the 
liquid flow dynamics can be expressed as a first order equation. This 
approach was used in this study, to express the liquid flow perturba- 
tions in the form: 


{ ' “t/t, 


Pee eee a (3.7-1) 


Tne T, is a liquid flow time constant and K accounts for different 
flowrate changes on each tray. 
Equation (3.7-1) can be expressed as the following differ- 


ential equation: 
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eas Wed (3.752) 


This expression relates the flow out of a stage to the flow into that 
Stage as a function of time. It must, however, be modified for the 
feed tray, to account for changes in feed flow. The equation for the 


feed tray is: 


De i ! 4 ‘ 
Le = 1, bt Ke 7 Kee = Lo] (3.7-3) 


where Kee accounts for different changes in F and Lp. 
The vector L also includes a perturbation term for changes in 
distillate flow since changes in this variable are observed for changes 
in feed and reflux flowrates. The response of distillate flow to feed 
flow changes can be written as: 
' ' -t/tp 
Oa ae [l-e ] (3.7-4) 
The response to reflux flow changes is slightly more complicated, since 
the distillate is withdrawn under condenser level control. A step 
change in reflux flow will produce an immediate equal and opposite step 
in distillate flow followed by a gradual response to a final steady 
State value. This can be expressed mathematically as: 
~t/t) 
D = “Lig + Lig Kpp [l -e J (3.7-5) 
If the tp terms in equations (3.7-4) and (3.7-5) are assumed to be 
approximately the same, the overall equation can be written as: 
-t/ -t/tp 
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This can be expressed as a differential equation in the form: 


e 4 Z be i] t - i] 1 2 

D = tp he KEp + Ligken = LioJ (3.7-6) 
Applying equation (3.7-2) to the appropriate stages and combining these 
with equations (3.7-3) and (3.7-6), the following matrix-vector differ- 


ential equation can be written: 
Pa) eb + ND Se) 


where the vector D includes feed and reflux flow disturbances, and M 


and N are appropriate coefficient matrices. 


3.8 Augmented State Space Model 
In the previous sections several expressions which can ve 
used to simplify equation (3.5-3) have been presented. Each of these 


equations is rewritten below: 


Wo + Wh = ACh + (Ah - BAIL + BUH + CC, - Ae, 7D (3.5-3) 
H = By h (3.6-1) 
L = ML+ND (eee!) 


Substitution of equation (3.6-1) into equation (3.5-3) and simplifica- 


tion yields: 


(3.8-1) 
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Equation (3.7-7) can now be combined with equation (3.8-1) to give the 


augmented state space equation: 


h ees |e ath ericc 
= - D (3.8-2) 
L geteH P.|Lt} ch 


where the submatrices are as defined previously and 0 is a 10 x 10 null 


matrix. This equation can be expressed in the form: 
x = Gx+tHD (3.8-3) 


and as such is a general state space representation of the distillation 
column. 

Equation (3.8-3) represents a matrix-vector differential 
equation of order 20. There are 20 states or time varying unknowns and 
5 inputs. These inputs appear in the vector D which comprises a 
control/disturbance vector. The model in the form of equation (3.8-3) 
can be solved for a specific D to yield the time domain response of 


each variable or state. 


3.9 Disturbance-Control Vector 
The possible disturbance and control variables for the dis- 
tillation column have been introduced into the model by the vector D. 


This vector contains the following perturbation inputs: 
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The disturbances which can enter by the feed stream are feed flow and 
feed enthalpy. Since the cooling water flow is used to control the 
nressure in the condenser, its effect as a load or disturbance vari- 
able must enter into the model and is included as Q.> a change in 
condenser heat duty. The two control yariables most often used, and 
used in this study are reflux flow and steam flow (reboiler duty Qn) 
and these enter into the vector D. This formulation allows the dis- 
turpance and control variables to be isolated from each other, and 


from the state variables. 
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CHAPTER IV 
MODEL EVALUATION 


4.1 Introduction 

Before a process model can be accepted as representative of 
a process, its ability to predict the effect of disturbances on the 
process must be verified. The simplest form of disturbance to utilize 
is a step disturbance since it is easily obtained experimentally and 
analytically. The approach adopted in this study is to obtain the 
time domain solution to the linearized state space model for a step 
disturbance. This time domain solution is then compared to the 
experimental response to the same step disturbance. This testing 
procedure, and the solution to the state equation are outlined in 


the following sections. 


4.2 Solution to the State Equation 
The general state space equation, in differential form, can 


be written as: 
x = GxtHD (4.2-1) 


where K 1s an nox |) state yector 

Gis ann x n matrix 

Dissan mx 1 disturbance-contro| Vectoy 

His ann x m matrix 
The solution to this matrix differential equation can be found in 
standard texts on state space techniques and has been outlined by 


Newel] [32]. The solution can be written as: 
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x(t) = (ty) x(ty) +f F(tyoe) HOC) de (4.2-2a) 
O 


where f(t) is a matrix function defined as 


F(a) = exp(G a) (4.2-2b) 


If the vector D is assumed to be constant over the interval {t; ty 


Te aes ths then equation (4.2-2a) can be written as: 


x(t;) = ¢ x(t.) + a D(t,) (4.2-3) 
where 
$ exp [G t,] 
a 
Ath | Beprenesse utd 
2 pala i 


This formulation gives an exact solution to equation (4.2-1) 
and can be used to yield the time domain solution. This approach can 
de programmed for a digital computer and a series approach used. However, 
if the model is of a large order the computational time can be quite 
large. An alternate approach whereby a trapezoidal approximation is 
used will be outlined below and the resulting solution compared to the 
exact solution obtained by solving equation (4.2-3). 

Consider a time interval of duration T and an incremental 
time At such that: 


T = k At 


If D can be considered constant over a time interval {t; oe EU ty ae 


then within this interval: 
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Also from equation (4.2-1) several general expressions can be written: 
Kit) eee Baxter HD(t) (4.2-4a) 
and x(t.) = ex(c) Het.) (4.2-4b) 


However in the interval {t; t) < t < t, + T}, D is constant and for 


2 < k and es 0, equation (4.2-4a) can be written as: 
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and for the following interval: 


x([e+1] t) = @ x([2+1]at) 
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Also if x(zAt) is known, x([z+1]at) can be calculated from the mean 


value theorem as 


XV IAt) = x(at) fat n (4.2-5) 
where ae) 
and gat Sa S$ [2+] Jat 


The trapezoidal approximation that: 
n = > Ck (Let Jat) + i(2at)] 


can pe made to simplify the solution. 
Application of this approximation and substitution into equation (4.2-5) 


yields: 


X( (s+) Jat) KCAL) + “(6 Migaty + ex(eel lat) + 2 HD] 


This equation can be rearranged to yield: 
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where I is ann x n identity matrix. This equation can be rewritten 


in the form: 


x(L2t1 Jat) = @ x(zat) + a D (4.2-6a) 

where 
e = f1-Ser' tr+S) (4..2-60) 
a= (1-Ahe6y' wat (4.2-6c) 


These equations can be used to obtain an approximate solution and the 
difference between this solution and the exact solution will decrease 
as At is decreased or as k is increased. It can be shown that, for a 
given input disturbance, the final steady state predicted by equations 
(4.2-6a) and (4.2-1) are identical. From equation (4.2-6a) the final 


value is obtained by setting: 
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= [I cs oe A D (4.2-7a) 
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From equation (4.2-1), 1f the system 1s stable, that 1s 17 the real 
parts of the eigenvalues of G are negative, the final steady state can 


be obtained by setting 


yielding PS ge g H D (4.2-7b) 
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The equivalence of these final values can be shown as follows: 
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This result is the same as in equation (4.2-7b) indicating that the final 


Steady state values are identical. 
The equations (4.2-6a, b, c) represent an approximate solution 


to equation (4.2-1). The exact solution is as given in equations 
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can be shown, and so some criterion for the choice of At given. 
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The similarities and differences between the two solutions 


The 


following is a presentation of this in scalar form for simplicity, and 


with the extension to the vector-matrix form. 


For the scalar case: 


the exact solution is 


hd dt (4.2-8) 


An approximate solution, by analogy to equations (4.2-6a, b, c) is given 


by 
x(t, ) ee Oa ) ea 
where 
ite t 
a= (1- = 9) ie = 9) 
t 
b = (l- = 9) \ ty h 


(4.2-9a) 


(4.2-9b) 


(4.2-9c) 


The time difference ty in this case is analogous to the At term defined 


for the matrix case. 


By the binomial expansion: 


So the terms in equation (4.2-9) can be redefined as: 
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; “9° co ec 
apa let Cea a at aCe (4.2-10a) 
ght,” g’nt,° gknt (k+1) 
ban tyh arte ies gs nes arr (4.2-10b) 


However, the coefficients as defined in equation (4.2-8) can also be 


expressed in series form as: 
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At this point, a criterion can be selected for the choice of t)> the time 
interval used. A convenient form, which expresses the error in the 
approximation as a fraction of the previous term in the expansion is 
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3.3 (ie 
t, g t g tig 
and 
ty“g*h ty “hg tig 
1? < B 5 or eas © 48 (4.2-11b) 


where a and 6 are selected constants. Since equations (4.2-lla, b) are 


essentially the same, a criterion can be selected in the form. 


a4 ae ie) (4.2-12) 


where a iS a measure of the goodness of the approximation. 
This criterion applies to the scalar case and the extension to 
the matrix-vector case is based on the definition of an eigenvalue. If 


the original equation (4.2-1) is diagonalized by the similarity transform: 


to yield 


then each canonical state can be expressed in the form of a scalar equa- 
tion. So for each canonical state a criterion of the form of equation 
(4.2-12) can be written. The logical choice of a At is therefore to 


select the smallest one as given by: 
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max [ te Re 04) ltaved. tan (422213) 


where the he are the eigenvalues of the matrix G. 
The time base of ¢ and A as in equation (4.2-3) is then at. 


The time base or sample interval can be changed to any time T if 


T= Kk ht 
by the relations: 
k 
ig, oe oe 
k-] . 
aiees oy A 
j=0 


to give the discrete form with time base T as: 


4.3 Experimental Studies 
4.3.] Experimental Equipment 

The distillation column used in this study is located in the 
Department of Chemical and Petroleum Engineering at the University of 
Alberta. The glass column is nine inches in diameter and has eight 
trays, each fitted with four 2-1/4 x 1-7/8 inch diameter bubble caps, 
1-1/2 inch diameter inlet and outlet wiers and a 1 inch diameter down- 
comer. The vapour is condensed in a glass total condenser fitted with 
stainless steel] cooling coils. The reboiler is a 7 inch nominal dia- 


meter basket type reboiler using 70 psi steam. 
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The column is instrumented with conventional local analog 
controllers as well as being interfaced to an IBM 1800 process control 
computer. The present configuration allows the computer to act in a 
Supervisory capacity by varying the setpoints of the local controllers 
for feed, steam and reflux flow rates. The computer also acts in a 
data gathering and logging capacity by sensing many other variables. 
Details of the column and its instrumentation have been adequately pre- 
sented by Svrcek [46] and by Berry [6]. 

The condenser outlet fluid composition is continually moni- 
tored by means of an in-line capacitance analyzer. The details of this 
apparatus have been presented by Svrcek [46]. The reboiler bottoms 
composition is monitored by the use of an on-line gas chromatograph 
operating under computer supervision. The separation and analysis is 
done on a 150 second cycle. The computer collects the data and analyzes 
it to calculate the composition. The values are also collected in a 
computer disk file for data logging purposes. 

The feed is pumped onto the fourth tray (from the bottom) 
under flowrate control. It is preheated by a steam preheater, under 
temperature contro]. The bottom product is withdrawn under reboiler 
level control while the top product is withdrawn under condenser level 
control. The reflux is preheated by a steam preheater before returning 
to the top tray, and it is maintained close to its boiling point. 

4.3-2 Open-Loop Jests 

Experimental transient data for use in the eyaluation of the 
proposed mode] was gathered by carrying out several ‘open-loop' step 
tests on the column. These tests were open-loop in that the effect of 


a step change in one variable on the product compositions was investigated 
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with other variables maintained at their steady state values. For 
example, if a step in feed flowrate was used, steam and reflux flow- 
rates, column pressure, reboiler and condenser levels, and feed and 
reflux temperatures were controlled at their steady state values. 

The step input disturbances were implemented by introducing 
a step change in the setpoint of the appropriate local controller by 
means of the IBM direct digital] control (DDC) system. The transient 
data was then gathered using the data accumulation algorithm in the 
DDC system. The data was stored in disk files and so could be plotted 


or punched onto cards for future use. 


4.4 Comparison of Simulated and Experimenta] Results 

The discrete form of the model, as represented by equation 
(4.2-6a) was used to calculate the transient model response to given 
inputs. The time base was taken as one minute for convenience. The 
disturbances, as calculated from the average of initial and final values 
were input to the model. In general, one of feed, steam and reflux flow 
was perturbed by a step in flowrate, and the others held constant. Since 
the cooling water flow was used to control column pressure, the condenser 
heat duty varied, and so was another disturbance input to the model. 

The overhead and bottom product enthalpies were calculated 
from the transient composition data. The bottom product analyzer, the 
process gas chromatograph, gave results with appreciable noise. An 
exponential filter was used to smooth the data. This, coupled with the 
measurement time lag, made a comparison of bottom enthalpies difficult. 

The following sections detail] the transients for various forms 
of step disturbances. Table 4.4-1 indicates the final values of the 


overnead enthalpy perturbations for the steps considered. Table 4.4-2 
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presents the same information for the bottoms enthalpy perturbations. 
In most cases the bottoms enthalpy transients are not presented, due 
to inaccurate measurements as described above. 

4.4-] Feed Flow Disturbances 

The ability of the model to adequately predict the effect of 
a feed flow disturbance on the experimental system was studied. The 
magnitude of the step changes used and the final steady state perturba- 
tions obtained are indicated in Tables 4.4-] and 4.4-2. The model and 
experimental responses of the overhead enthalpy for feed changes are 
presented in Figures 4.4-] and 4.4-2. A similar set of data for the 
bottoms enthalpy is presented in Figure 4.4-3. These data are presented 
as a fraction of the total change complete, relative to the final model 
steady state. 

In the case of the positive feed flow step, the tendency is 
to increase the overhead composition. This tends to force the composi- 
tion into the equilibrium curve pinch area. The change observed is 
therefore smaller than that predicted by the linearized model. The 
model responds slightly faster also, and a delay is observed between the 
time of disturbance and the experimental response. Figure 4.4-3 shows 
the bottoms enthalpy response to the same feed flow disturbance. The 
experimental points illustrate the amount of measurement noise in the 
filtered signal from the gas chromatograph system. The final steady 
states, as indicated in Table 4.4-2 show moderate agreement, but the 
transients do not agree very well. 

The overhead enthalpy response to a negative feed flow step 
are presented in Figure 4.4-2. The model predicts a smaller change 


than experimentally observed, as the composition moved away from the 
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equilibrium curve pinch area. It can be seen that, again, the experi- 
mental curve exhibits a short time lag. 
4.4-2 Reflux Flow Disturbances 

The ability of the model to predict the open loop response to 
reflux flow changes was tested next. The magnitude of the step changes 
are presented in Table 4.4-1. Figures 4.4-4 and 4.4-5 indicate the 
transients obtained in each case. The figures indicate the overhead 
product enthalpy as a fraction of the final model steady state value. 

In the case of a positive flow change, the tendency is toward 
an increase in product composition. This drives the composition toward 
the pinch area of the equilibrium curve so the model predicts a larger 
perturbation than was actually observed. The model also exhibits an 
initial] time lag behind the experimental curve, but it is not very large. 

For the case of a negative step in reflux flow as indicated in 
Figure 4.4-5, the agreement is quite good. The final value predicted by 
the model is slightly below the experimental one, as the composition 
moved away from the pinch area. Again the model responds slightly slower 
initially, so the model response lies entirely below the experimental 
one. 

The final yalues of the bottoms enthalpy perturbations for the 
two cases above are presented in Table 4.4-2. 

4.4-3 Steam Flow Disturbances 

The experimental response of the overhead enthalpy perturba- 
tions to steam flow changes are indicated in Figures 4.4-6 and 4.4-7. 

The responses are presented as a fraction complete, relative to the final 
model value. The step in steam flow caused a step in reboiler duty which 


in each case was accompanied by a change in condenser duty of the same 


- 
ef 
mn 


vant node ah neo St ab Han 
ed "palnaes Las hee & 20 td a 
bw | a ee 
rsancewse rh WOE oP A = P 


otps1g. of fabout 2° 1D uth. 


at Sanogesi qool neqe a) 


' t ; & =)" =a Ath 7 ; E : 
zopnetio gate O17 to Sbudiaos/ ont .2xSn botest caw 2SpNSiD we 


yh +hdeT nf beaded 


. " : i 4 witiret | s. 
wis SoeoTbAD C-h.P bite o-hee e2TUR I)! 


=e 
0, 
Seotevo sdt-aieotuul. esvupi? S(T .senz foss MU borirsi do asost aad 
7 Mt 
; _ ro on. Ponts + 9 idoBNt Be 2 wilt Jouborg 
»SUTEY SR BF 2 Yok Labom Bt SI FO Wits pte 2k ug! f 8 7 


biawod er yousbret of) ,Sbasrs wort syiireog & 1%, g26> sit ni > a 
a 
hyswot ftottieognos adi zeylal 2tt) nots r2oqmos Jounorg vi 12698 a5 


; "spt6! & etorpony [ob it Of SVD mii?! Tepe ore To ss¥p ante 


> 
, ‘ ~~ man 5 rep cif ae ai rd 
als eT radrexs Heise isp “it ] ba 19200 Vii bUl gS caw ned ; norte aad 


ees 
«3p iBT y 19V FOU 2 cs + [ 7] . lg he) 1 oe ne ah . 7 fia hnfaed pel ne ae late me 


at batexibtir ‘26 wolt *ultsy oi uste avitapon 6.70 Bebo Bn TOF 


fant? Sai oop attup =f somes pr S01 -e-h aeS 


Ke Phat : " y 
yo GOJDTbSIG Siti py , 


nore feoum oD ant ely P TATYE) bz heel ait) 2m ah wofsd Vise; fz ef | ab va Pa 


i 


“auote Visirotfe ebnouesy fabum sit-nieg! .ssys dontq, sat meqy oe beve 


fetnemiisaxe ond woled vieitits est| senogesy [shom 5m0 ee vais oa 


( 
i-h.b oldsl wi bedaseore e165, 9v0ds gedeo OM 
~~ a 7 


t-b : 
- 


Sy (Ot. shore acyu3 sq vats Jus metiod- srt ta e5uley Tent? SAT. mt 
oa 


sansdyur e til word npore: 


-pdywwiya9g. vgrtsitine beofisvo sit to senoqes (sinonivaqsa’ oat 7 


.\-b.8 bos a-Bph 2quup it ar betasibar 916 eopnena "0 aos: or 
r 7 7) : 7 


teatt arid of avia6! af SFal gino’ aves fe eT 25 20g 
ba tein tub ast igne, ni qa? sé bozue: wort “M63 92 2 i qote = = a 


. 


eee t, 
- waciead eos 5 vaeiiabno. i ape sii BNO by 
ica ak 


QVERHEAD ENTHALPY 


REFLEX FLOW 


RS A Seaipaiense?: Wine Were Berane = SOR Cee eth ont coer eco ee ee OME Vein ede OPE BD ALY 1 ESP 6 Sv eee 


FIGURE 4.4-4: 


20 40 GO BO 
TOE PN APTS 


COMPARISON OF MODEL AND EXPERIMENTAL OVERHEAD ENTHALPY 
RESPONSE TO A POSITIVE REFLUX FLOW STEP CHANGE 


400 


wh st ie 


m Hg 13, gtaVvO. iAT¥ON via A = ioe 
“omni 1) Saco fini ou on ee 


OVERHEAD 


S + EXPER, eerie eae 
ua == Menen ve: 
CL + 
= 2 
as + 
Pa 
Lil . 

+ 
ate 
+ 
+ 
OQ 


REPELUX fLeW 


AQ (50 G0 400 
Tai” TING MIND es 


FIGURE 4.4-5: COMPARISON OF MODEL AND EXPERIMENTAL OVERHEAD ENTHALPY 
RESPONSE TO A NEGATIVE REFLUX FLOW STEP CHANGE 


al 


a Mit Sar 
nia 
Jey 
a 


ns 
TH ; vr 

i Take VO. at 

— aK 

19. AKT) 304 ! - 

Aen 


ces 


OQVERHEADQ ENTHA 


STEAM FLOW 


yy 


- + EXPER: 
— MODEL by oct a a oO 
“hE 
+ 
+. 
we 
+ 
+ 
“+ 
ob 
(@) 
Q 
(u 
fo | 
BO La al Se ben 
O 
DX 


d 


a 
O 20 AD GO BO TOO 


Pie UN MINGES 


FIGURE 4.4-6: COMPARISON OF MODEL AND EXPERIMENTAL OVERHEAD ENTHALPY 


RESPONSE TO A POSITIVE STEAM FLOW STEP CHANGE 


; j=: See i. 
ree roe. fn pel — ae Le eprom 


~~ 


ar oe : i a — 
_ ry aaa 7 


ator 4 


 § SATU ON: ik 1703 ns A 2° 


vw: A ¢ 1 4) 
| 14 = sail te 10 ae i 


BVERREAL  ENTAALPY 


STEAM FLOW 


FIGURE 4.4-7:; 


Tv IN IV SIS 


COMPARISON OF MODEL AND EXPERIMENTAL OVERHEAD ENTHALPY 
RESPONSE TO A NEGATIVE STEAM FLOW STEP CHANGE 


54 


order of magnitude. 

For both the positive and negative step changes, the model 
response compares very well with the experimental one. The model 
reaches a final value quite rapidly, probably due to the assumption of 
no vapour flow changes. 

Again, the final yalues of the bottoms enthalpy perturbations 
for these step inputs are presented in Table 4.4-2, but the transients 


are not shown. 
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55 
CHAPTER V: MODEL REDUCTION 


5.1 Introduction 

The model of the distillation system under study as formulated 
in Chapter III of this work is of order 20. In general, the order of 
such a model will be 2(N+2), where N is the number of trays in the 
distillation column. For a control study on a large system, the Sie 
of the problem thus becomes prohibitive. Also the usual control 
strategies require values for all states, and in general, many of the 
states in this model cannot be directly measured. So the feasibility 
of model reduction was considered. Two approaches were considered and 
their results compared to the original models and to each other, as 


outlined in the following sections. 


5.2 Marshall's Method 
5.2-] Reduction Theory 

This method is the so-called "modal approach" outlined by 
Marshall [31]. The method has been further studied by Wilson et al 
[53] and extended to the case of a discrete model. The continuous 
method, as implemented in the GEMSCOPE system [54] was used in this 
Study. 

The method revolves around the partitioning of the states 
into fast and slow mode states. The reduced order model is comprised 
of the slow states, and an instantaneous response is assumed for the 
fast canonical states. The method reduces a model of order n of the 


form: 


y¥ = Ax +Bu (5.2-1) 
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56 
to the reduced model: 
4 7 Bp hy FBR (ee 


where Xx, are the states to be 'retained', and Ap and Bp are derived from 
A and B as will be indicated below. 


Equation (5.2-1) can be rewritten in a partitioned form: 


= + u (5.2-1a) 
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where X, are the states to be retained, and X5 are those to be eliminated. 
An n xX n matrix M exists which will transform equation (5.2-1) into the 


Jordan canonical form [53] by using the similarity transform: 
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where z is the canonical state vector. If the states x, are to be the 
‘slow mode' states, matrix M is arranged so that its columns are ordered 
from left to right, in order of decreasing significance of the corres- 
ponding eigenvalue of the. matrix A. When equation (5.2-3) is introduced 


into equation (5.2-la), the Jordan canonical form results: 
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where a which is a block diagonal matrix, and 8 are defined by: 
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can be obtained from equation (5.2-3) and (5.2-5c), 


Soy = p, py IPs x. ] (5.2-7) 


Since the Z, canonical states are assumed to correspond to the fast 


modes, assume that: 
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This expression is then substituted into equation (5.2-7) to yield: 
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= -] -] 
Xo = Ba [- a Bo U- pz xX] 


This can then be substituted into equation (5.2-6) to yield: 
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For the extension of the above method to the discrete model case, and 
for equivalent forms of the Ap matrix see Wilson et al [53]. 
5.2-2 Reduction Results 

Marshall's reduction method as implemented in the GEMSCOPE 
system [54] was used in an attempt to reduce the 20 state distillation 
model, to one which could be used for a control study. This reduction 
method permits the choice of which states and which eigenvalues to 
retain in the reduced model. Several combinations of states and eigen- 
values were tried, in an attempt to obtain a valid reduced model. The 
problem was assumed to be one of removing all unmeasurable states, i.e. 
internal liquid flows, and as many enthalpies as possible. The models 
obtained were compared to the original one for a given disturbance and 
thus accepted or rejected. No acceptable two state model was found. 

One of the better models obtained was a fourth order model in 
which the states were the enthalpies in the first, fourth, sixth and 
tenth stages. Figures 5.2-1 to 5.2-3 compare the end product enthalpy 
perturbation transients in normalized form as predicted by the four 
and twenty state models, for various step inputs. 


In Figure 5.2-1 the responses to a feed flow disturbance are 
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shown, and good agreement between the models is observed. In the case 
of a reflux flow step, illustrated in Figure 5.2-2 the four state model 
exhibits a response that is appreciably faster than the original large 
model. For the reboiler duty disturbance, illustrated in Figure: 552-3, 
both models exhibit rapid responses, with the four state model 
exhibiting some overshoot. In all cases, the final steady state values 
are comparable since the reduction method preserves the final steady 
State values, for step input disturbances. 

As stated earlier, several combinations of states and eigen- 
values were tried, in an attempt to obtain an adequate model. Many 
problems were encountered in these attempts to reduce the model, such 
‘aS "instant" responses and "infinite" final steady state values. It 
is thought that these problems stem from the grouping of eigenvalues, 
and from selecting some from each group. Table 5.2-1 presents a tabu- 
lation of a sample set of eigenvalues of the 20 state model. 

Because of these problems, and the desirability of a two 
state model, in which only the product enthalpies appear, it was decided 
to attempt another, more analytical reduction method. This method is 


outlined in the following sections. 


5.3 Anderson's Method 
5.3-1 Reduction Method 

This method is a ‘least-squares approach' to model reduction. 
The method has been outlined by Anderson [2] and is used to find a 
reduced discrete model which approximates the response of a larger model. 
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a reduced model: 


x, (k+1) = op x, (k) + Ap u(k) (5.3-2) 


can be found, which approximates the response of the original model for 
the states to be retained, designated as x,» the first r states. 

The method requires the values for the states and disturbances 
for a number of time intervals, say (2+1), where 2 is the number of 


points to be fit. Form the vector: 


by = [%q (1) E xg(2T) bee fxg (G41) 


where T is the discrete model time interval, and q refers to a specific 
state. Also the vector om can be formed of the elements of dp and dp 


to be determined as: 


If the matrices B and 
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where the first r states are to be retained. 


Also, definition of a set of vectors of the form: 
Lr [x, (iT): xX» (iT) cae x,.(iT) u, (iT): fauna.) 


allows the formation of a matrix m such that: 
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By inspection of matrices B, 


HWO> 
m 
= 
[@E 


m it can be seen that the time domain 


solution of the first r states is: 


B=mc (5. 3-3) 
and that C will yield dp and Ap Since: 
C = [ep : ag] (5. 3-4) 


Anderson [2] has indicated that the solution does exist and also mini- 
mizes a quadratic error index of the least squares type. 
The solution to Equation (5.3-3) cannot be directly obtained 


and a pseudo-inverse approach is needed. This approach yields: 


WS > 


C! = Top: ap] = 8 


5.3-2 Reduction Results 

| The reduction method outlined in the previous section is 
essentially one of fitting a set model form to transient data generated 
for a set disturbance by a least squares technique. The simplest form 
of disturbance to use is a step disturbance, so this form of forcing 
was used. However, only one step could be used at a time or the matrix 
im! m] would be singular. Because of this, each disturbance was fit 


in turn, so five reduced state transition matrices ($p) and the five 
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vectors Sp.» the columns of Ap were obtained. Since the distiltation 
column responds differently for various disturbances, the five state 
transition matrices obtained are different. 

The matrices corresponding to the feed flow, reflux flow and 
feed enthalpy disturbances were all quite similar. The matrices 
corresponding to reboiler and condenser duty were similar to one another 
but were quite different from the other three. Since one matrix dp 
must be used, one very similar to the former three was chosen, and the 
columns of the matrix Ap manipulated to give the same final steady 
States as the original reduced models. The effect of this approximation 
is to slightly slow down the responses to condenser and to reboiler duty 
changes, while keeping the responses to the other three disturbances 
essentially the same as the original model. The resulting reduced model 
is presented in Figure 5.3-1, with a sample interval of one minute. 

The responses of this two state model to several disturbances 
are compared to those of the twenty state model in Figures 5.3-2 to 
5.3-4. The responses to a feed flow step are very much alike, as 
indicated in Figure 5.3-2. The responses to a reflux flow step are 
shown in Figure 5.3-3 and it can be seen that the reduced model responds 
faster than the original model and does not exhibit an initial lag. 
Figure 5.3-4 compares the responses to reboiler duty changes and it can 
be seen that the reduced model is appreciably slower than the original 
model. The application of this reduced model, therefore, has the 


effect of slowing down the response to a steam flow change. 


5.4 Conclusions 


The reduced models described in the previous sections both 
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compare very well to the original model for given disturbances. The 
four state model compares better than the two state one for reboiler 
duty changes , but the converse is true for reflux flow changes. The 
two state model has the obvious advantage that both states can be 
easily calculated from the measured end product compositions, whereas 
the four state model has two states, the fourth and sixth stage 
enthalpies which can only be approximated from noisy temperature 
measurements. Because of this disadvantage and the comparable response 
of each model, it was decided to use the two state model for control 
studies. 

Before proceeding to control studies, the two state model 
was compared to the experimental data used to test the original model. 
Figures 5.4-1 to 5.4-4 show the comparison of the model and the experi- 
mental data. In each case the model response is comparable to the 
original 20 state model as presented in Chapter IV, and the comparison 


to experimental data shows moderate agreement. 
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CHAPTER VI 
CONTROL STUDIES 


6.1 Introduction 

The objective of this study was to attempt to formulate a 
control strategy whereby both the overhead and bottom product com- 
positions would be simultaneously controlled at preset values. Ina 
Single loop configuration the reflux flow is used to control top 
product composition and heat input to the reboiler, or steam flow is 
used to control the bottom product composition. However, as described 
previously,the reflux flow affects the bottoms composition, and the 
steam flow affects the top product composition. This interaction can 
give rise to problems when simultaneous control of both product com- 
positions is attempted. 

In any attempt to contro] the end compositions, some indica- 
tion of their actual values and variations must be available. In this 
Study the compositions themselves were measured and used to adjust the 
specific control variables. The overhead product composition was con- 
tinually monitored by an in-line capacitance cell, and the bottom com- 
position by a gas chromatograph operating under digital computer 
Supervision. 

In the control studies, tle first step inyolved performing 
several single loop tests, changing the feedback gains, to establish 
"Suitable" values. This was done for both top and bottom product 
control loops. The gains found in these tests were then combined in a 
multiloop configuration to demonstrate the extent of interaction 
between the loops. The state space model was then used to derive a 


multivariable control scheme for experimental evaluation. 
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6.2 Single Loop Control 

As previously mentioned the first step in the series of con- 
trol studies was to perform single loop tests. This served two pur- 
poses; to determine actual feedback controller gains, and to test the 
digital computer control programs. However, prior to the experimental 
tests, the controllers were studied by digital simulation using the 
two state model. The following sections describe these tests and 
their results. 
6.2-] Simulation Studies 


The model and general feedback control Jaw can be designated: 


x(k+1) = 4 x(k) + 4 d(k) + 8 uk) 


For the single loop tests the feedback matrix Kep contained only one 
element, the one completing a single composition control loop. Figures 
6.2-] to 6.2-3 show the open loop, top composition controlled and 
bottom composition controlled simulation runs respectively. The top 
feedback gain used was 22 and the bottom gain was -330. The simulation 
results show that these yalues nearly result in unstable operation, so 
they should be the upper limits of possible yalues. As the gains are 
decreased (absolutely), the offsets increase and oscillations decrease. 
It can be seen from Figure 6.2-2 that the top composition control loop 
has a large effect on the bottom composition, and actually has a con- 
trolling effect since it decreases the offset. 

In preliminary experimental tests, it was noted that there 


was measurement noise in the bottoms composition from the gas 
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FIGURE 6.23: SIMULATED BOTTOM ENTHALPY CONTROL BEHAVIOR (Kop = -330) 
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chromatograph system. It was decided to introduce this noise into the 
Simulation and to add an exponential filter to minimize its effect. 

The resulting signal was obtained by superimposing an experimental ly 
measured noise signal which had been filtered onto the simulated signal. 


The filter used was an exponential filter which has the form: 


Xe(k+1) = 0 x(k) + (1 - a) x 6ktT) 


where x = actual measurement 
X = filtered measurement 
a = filter constant 


Figures 6.2-4 to 6.2-6 present the simulation results corresponding to 
Figures 6.2-1 to 6.2-3 respectively, with the filtered noise added. 
The filter constant used was 0.8. The open loop gains had to be 
decreased since the noise tended to yield an oscillatory system. The 
bottom gain used was -300 and the top gain was 20, and in each case 
instability was approached. 
6.2-2 Top Product Composition Control 

The feedback controller gain as found by simulation was then 
tested experimentally. The disturbance used was a feed flow step. 
The simulation gain of 20 was found to result in oscillatory operation. 
This is probably due to process noise. A value of 10 was then tried 
and found to give a stable result, and moderately good control. The 
results for both a positive and negative feed flow step are given in 
Figures 6.2-7 and 6.2-8 respectively. In each case there are only small 
fluctuations in top product composition. However, there is a very 
small offset in final steady state values as is expected for feedback 


control when a step input is used. 
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FIGURE 6.2-4: SIMULATED OPEN LOOP BEHAVIOR - FILTERED NOISE INCLUDED 
IN BOTTOM COMPOSITION MEASUREMENT (a = 0.8) 
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FIGURE 6.2-5: SIMULATED OVERHEAD ENTHALPY CONTROL BEHAVIOR - MEASUREMENT 
NOISE INCLUDED (q = 0.8, Koy = 20) 
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FIGURE 6.2-6: SIMULATED BOTTOM ENTHALPY CONTROL BEHAVIOR - MEASUREMENT 
NOISE INCLUDED (a = 0.8, Kop = -300) 
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The bottom composition response is also shown in these 
figures. In the case of Figure 6.2-7, the bottom composition reaches 
a final steady state offset of about 0.4%. This can be contrasted 
with an expected offset of 0.8 to 1.0% for a corresponding open loop 
case. This is the same controlling effect as predicted by the simu- 
lation results and illustrated in Figure 6.2-2. In the case of Figure 
6.2-8 the bottoms composition reaches a composition less than about 
0.4% where the bottoms gas chromatograph cannot measure it. 

6.2-3 Bottom Product Composition Control 

As outlined previously, there was considerable noise in the 
bottom composition measurement and so an exponential filter was added 
to the simulation results. A corresponding filter was added to the 
experimental control system using the same constant as was used in the 
Simulation tests. The control loop gain in the feedback matrix was 
adjusted to give a stable response, and moderately good control. The 
Simulation value of -300 was found to yield an unstable system, and 
so smaller values were tried. A gain of -50 was found to yield stable 
operation, whereas larger values tended to drive the system unstable if 
a large noise level appeared in the bottoms composition analysis. 

Figures 6.2-9 and 6.2-10 present experimental results using 
this gain, for a positive and a negative feed flow disturbance 
respectively. It can be seen from Figure 6.2-9 that the bottoms com- 
position exhibits appreciable noise and an offset of about 0.3% com- 
pared to the comparable open loop test where an offset of 0.8 to 1.0% 
is expected. This is a definite improvement, however, the offset must 
be tolerated due to the necessity of using a low gain value. In the 


case of Figure 6.2-10, a similar but opposite result is observed with 
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a resulting offset in bottom composition of about 0.5%. In both cases 
the top product composition varies due to changes in steam flow. This 


effect can be seen most markedly in Figure 6.2-10. 


6.3 Multiloop Control 

The extension of the previous case, that is two single control 
loops, is to operate both loops simultaneously. This yields a two point 
composition control system. This approach was tested by both simulation 
and experimental testing, and the following sections outline the results 
Of eacn test. 
6.3-1 Simulation Results 

The two point control system was simulated using the two state 
model with measurement noise and filtering added. Initially, the con- 
troller gains as found in the single loop simulation studies were applied 
and were found to result in unstable operation. These gains of {-300, 
+20} were then tuned, and values of {-150, +15} were found to give a 
stable system. Figure 6.3-1 presents the controlled responses for these 
controller settings, The interaction between the control loops can be 
seen by the similarity in the shapes of tne reflux and steam flow curves. 
6.3-2 Experimental Results 

The two point control system was then tested experimentally. 
Tne composition feedback controller gains from the single loop tests were 
{-50, +10} for the bottom and top composition control loops respectively. 
These gains were used in the two point system, and feed flow steps were 
used to disturb the column. Figures 6.3-2 and 6.3-3 present the con- 
trolled responses obtained. In Figure 6.3-2, the bottom composition 
shows large oscillations and is tending to be unstable. The overnead 


composition, however, shows smal] changes only, and appears wel] 
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FIGURE 6.3-2: EXPERIMENTAL MULTILOOP CONTROL BEHAVIOR: +5% FEED FLOW 
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controlled. Figure 6.3-3 shows the response to a negative feed flow 
step, and it is noted that the bottom composition is stable with large 
oscillations and an offset of about 0.7%. The overhead composition 
again shows small fluctuations but also has an offset which is larger 
than in the comparable single loop case. 

The difference in observed stability of these cases may be 

due to the magnitude of the bottoms composition change expected for 

the input disturbance. The change for a positive feed step is expected 
to be larger than that for the negative feed step. The large change in 
bottoms composition as a result of the positive feed flow step could 
result in a shortlived period of control yielding instability. 

From these multiloop tests the following conclusions can be 

drawn: 

i) The controller loop gains in a single loop system may 
yield an unstable control system when applied simul tan- 
eous ly. 

ii) The control loops chosen for this study exhibit inter- 
active effects when applied simultaneously. 
iii) The interaction can be predicted by the two state model 


of the distillation column. 


6.4 Multivariable Control] 
6.4-1 Multivariable and Multiloop Control 

In a multiloop control system there are several control loops 
operating independently of each other. However, as demonstrated in the 
previous section, the control actions of these loops can interact and 
lead to a deterioration in control. The extension to a multivariable 


contro] system allows the loops to interact and the control action taken 
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compensates for this interaction. This approach should therefore lead 
to better control performance than the simple multiloop case. 
6.4-2 Multivariable Control Laws 

The application of a multivariable control system requires 
the availability of a process model. The most widely used form is the 


stationary state space form: 
X= Wee Did Bou (6.4-1) 


where n x 1 state vector 


~< 

iT] 
m 
= 


m x 1 disturbance vector 


oO 
i} 
re) 
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2 x | control vector 
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A, B, C = time invariant coefficient matrices of appropriate 
dimensions 


This model can also be expressed in the discrete form: 


x(k + 1) = $ x(k) + 4 d(k) + B u(k) (6.4-2) 


where k refers to a specific time interval and the matrices and vectors 
are the discrete forms of those in Equation (6.4-1). 

Many types of control laws have been proposed in the litera- 
ture. The best known are so-called optional laws derived for some 
specific form of performance criterion. Some of these optional laws 
have been presented by Newell [32] and by Hu [22], both of whom have 
applied a dynamic programming approach. In this study it was decided 
to first employ simple suboptimal control, followed possibly at a later 
stage in the study by optimal control. 


A simple suboptimal law can be derived from the discrete model: 


x(k + 1) = o x(k) + A d(k) + 8 u(k) (6.4-2) 


_ 


_— y wat sree veal iii ‘abit neon ta itt 10 > 


(T-h.d) et bU+ 2 
WILY ; te. = A 
os pay ! bead nm = ob 
e 
"“Qio5v + Oudno x mh = v 
: : “—_ 7-29 ah Pe T t Ss \ | 
aiarrdoims 10 220Pry16m Jet f— THB PTE » A 
' : a 
clld! ensamnt a — 


[ane A _ 7 
a 7 ; ay ~ i = an 


* 
é 


~2460- goo! itum stiméz silt esis sooknoaa: Fonsi * vedas 
_ = : 
pws! forhnod staaiaen TUM S-b. 


. eprtupey mateve fortaos ofdatisvidtun b 46 notesotldgs ait 


7 
odd 2) pict boo. yiebtwaecom art «| sbom esc 6 to yoil tdettev sw 


“wot somje Shate« 


_ 
= - 
ot aysTv2eib slg. a! bezzeyqes Sd gels ime fobon 2hat- 


(Sah. 9) ‘ at) F }. ( 1)b A (aye is = if ih WK - 
- 7 7 
— b 7 
Syotosy bas! 2aotasnW adi DAS Tevindil amit ahvissqe © of e7eree A eee 
(1 *) 0 AOPT Gall | Semi) Fo ny tor sjeroerb dt hs , one 


-pisi(f ott af Hbezono, SAG OVE H-.2Wwel !Tossrio> 24 eogys yeh i : 


auntOe 407 bav mob wel. fenptito bafiano-oe Sek rein o5t. af wt 

4 4 7 —_—" : uv oy 7 PE ane 

ewil Tenoitqa Seat to omod ~ .nurvstina fonemotyaq 70 :ane “a 999 

avail monw to tod , bes] uh Vd bne-[ SE) Tfswet a a 93 

bsbiosb esw tt gous ai Al .ARCIQAHt or Home hice > 6 bafT ge 

qatel 6 tn yldbezaq bpotton air ve (smtaqodua sgn olan tan 
_ ; a 

92 ieayeae lamrrqo xd Ybuse si st pet; 


se AS dere Boe 
9d 6 bY PAE 


piaatb it mot 


pA 


Wagers 


.. 
_ * 


96 


The objective of the control action u(k) is to drive the values of the 
states at the next interval, KC KET) to the state setpoints for any 
value of x(k) and d(k). Let the setpoint of the state vector by Xo 
and assume that the controls are unconstrained, that is, the required 
changes in controls are physically realizable. Then the control action 


will make: 


SAREE axe 


The required control action can then be found by solving for u(k) in 


the expression: 
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Since, in general, 8 is not square (n # 2), a pseudo-inverse approach 
[19] must be used. The control law can be expressed as a feedback- 


forward-setpoint control law: 


=sp 7s 
where Kea = -[g'g] | a’ 
Ses = -[e'e gia 
Ke» = Oster!" 


Figure 6.4-1 presents the two state discrete model of the distillation 
column, and the resulting control law derived as indicated above. 
6.4-3 Simulation of Feedback Control Behavior 

The feedback matrix as defined above was used for simulation 
studies. The setpoint was assumed to be the origin in state space. The 


effect of neglecting the feedforward term is to cause an offset in 
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FIGURE 6.4-1: MODEL AND CONTROL LAW 
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output for a step input. This can be seen in Figure 6.4-2 where the 
ideal case is considered, i.e. no noise in measurements. Figure 6.4-3 
is a similar result with the noise and filter added where it can be 
seen that controlled responses show small offsets and oscillations. 

This feedback matrix was then applied to the experimental 
system, and was found to yield an unstable system. The reason for the 
instability was investigated and found to be partially the feedback 
element which adjusted the reflux flow for bottoms enthalpy changes. 

The noise in the bottoms composition measurement introduced oscillations 
in reflux flow, which made the top control loop unstable. Also the 
large gain in the bottom control loop caused large oscillations in 

steam flow. It was decided to reduce these two matrix elements. The 
bottom loop gain was reduced to -50, the value used in the multiloop 
tests. The other gain term, an interaction term was set to zero. This 
has the effect of making the reflux flow independent of bottom enthalpy. 
However, steam flow changes are still dependent on the overhead enthalpy 
changes. This removes part of the multivariable effect of the control 
system, however, as demonstrated earlier the reflux flow has an inherent 
tendency to control the bottom composition. Because of this, it was 
hoped that the above matrix adjustments would yield a behaviour only 
Slightly different from the ideal controlled response. However, this 
was not the case as can be seen by the results presented in Figure 6.4-4 
which shows the controlled response with the feedback matrix adjusted 

as described above. By comparing this figure with Figure 6.4-3 it can 
be seen that these modifications cause larger offsets in both composi- 
tions and larger variations in bottoms composition. 


As well as the large feedback gains in the original system, 
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another cause of the instability could be the time delays in the control 
system. The gas chromatograph separation and analysis cycle is an 
inherent measurement time delay in the system. The filtering of the 
bottoms composition measurement also introduces a form of time delay. 
Investigations showed that with the exponential filter used in this 
study, and a filter constant of 0.8, a maximum delay of about 5 sample 
intervals could be expected. Also when the filter constant is reduced 
to 0.5, the time delay is reduced to less than 2 sample intervals. 

This time delay was added to the simulation work, and some tests done 
to evaluate the filter constant which should be used experimentally. 

Figures 6.4-5 to 6.4-7 present the results of simulation tests 
showing the manner in which the value of the filter constant that is 
employed affects the control behavior. The effect of the total 
time delay has also been shown. Table 6.4-1 presents the filter con- 
stant, total time delay, and feedback matrix for these tests, as well 
as for those presented in Figures 6.4-2 to 6.4-4. The total time delay 
is less than 4 minutes when the filter constant is 0.5 and is about 6 
minutes when the filter constant is 0.8. 

Figure 6.4-5 is the simulation result corresponding to the 
experimental test which was found to be unstable. It can be seen that 
there are large sustained oscillations in bottoms composition, and that 
if a measurement was encountered which represented a noise level larger 
than that used in the simulation, it can be seen that the system would 
tend to be unstable. Figure 6.4-6 presents the simulation results with 
no time delay and a filter constant of 0.5. The effect of reducing the 
filter constant can be seen by comparing the response with that in 


Figure 6.4-4. This has the effect of increasing the noise level in the 
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FIGURE 6.4-5: SIMULATED MULTIVARIABLE FEEDBACK CONTROL BEHAVIOR WITH 
TIME DELAY AND MEASUREMENT NOISE INCLUDED (a = 0.8, 
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in the bottoms composition measurement, but also reduces the inherent 
time delay due to filtering. Figure 6.4-7 shows the simulated control 
behavior using a filter constant of 0.5 and a time delay of 4 minutes, 
which is about the maximum that would be expected experimentally. 

This transient can be compared with that in Figure 6.4-5 and it can be 
seen that the magnitude of the oscillations is reduced although the 
frequency is increased. This test is a good estimate of the conditions 
which should result experimentally and show a stable system with sus- 
tained oscillations. 

6.4-4 Experimental Feedback Control Behavior 

The simulation results indicated that a filter constant of 
0.5 should be used with the adjusted feedback matrix. This test, as 
indicated in Figure 6.4-7, was attempted experimentally. Figure 6.4-8 
presents the results of a test with no input disturbance, which was 
done to see if the noise level tended to yield instability. As can be 
seen, the system is stable, with moderate noise in the bottoms com- 
position. 

Figures 6.4-9 and 6.4-10 compare the controlled and uncon- 
trolled responses for a positive and a negative feed flow step respect- 
ively. As can be seen in Figure 6.4-9, the sustained oscillations as 
predicted by simulation are present. Also the controlled offset in 
bottoms composition is only about one-quarter of the uncontrolled off- 
set. The distillate composition shows oscillations in the controlled 
case, but almost no average offset. In the case of the results using 
a negative feed step, shown in Figure 6.4-10, the oscillations don't 
appear to be as drastic as in Figure 6.4-9. The controlled offset in 


bottoms composition, in this case is only about one-fifth of that in the 


uncontrolled test. 
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FIGURE 6.4-8: EXPERIMENTAL MULTIVARIABLE FEEDBACK CONTROL BEHAVIOR - 
NO DISTURBANCE (a = 0.5) 
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These experimental tests support the results predicted by the 
Simulations, that is a reduction in average offset of bottoms composi- 
tion with sustained oscillations. The overhead composition also shows 
small oscillations in response to the bottom loop control action, 
namely, steam flow rate changes. In general, the changes in composition 
are small, and indicate a satisfactory feedback control law, despite 
bottoms composition noise. 

6.4-5 Simulation of Feedback-Feedforward Control Behavior 
The feedback-feedforward control law as defined earlier has 


the form: 


When this control law, with the original matrices is applied to the 
model, with no noise, the result is perfect control, that is, no changes 
in composition for any disturbance. However, as indicated previously, 
the experimental system has both a measurement time delay and measure- 
ment noise. Also the feedback matrix was modified because of the 


bottoms composition noise to yield the modified control law where: 


0.00 1229 -1.19 -.0169 .0369 
=FB ape aa “FF ie -.0831 a 

This control law was applied to the model with noise and time 
delay terms added. Figure 6.4-11 shows the simulated controlled res- 
ponses to a positive step in feed flow. The feedforward action causes 
an initial change in control action which makes the composition change. 


There are several large spikes in bottoms composition initially which 


tend to damp out with time indicating that the feedforward gain may be 
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FIGURE 6.4-11: SIMULATED MULTIVARIABLE FEEDBACK-FEEDFORWARD CONTROL 
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too large. This hypothesis is Supported by comparing Figures 6.4-11 and 
6.4-7 in that an initial deterioration in control appears on inclusion 
of feedforward compensation. As in the case of feedback control only, 
there are sustained oxcillations in bottoms composition but the sys tem 
appears to be stable. The bottoms composition change has an average 
value of about zero, as expected. 
6.4-6 Experimental Feedback-Feedforward Control Behavior 

The feedback-feedforward control law as presented in the 
previous section was tested experimentally. Figures 6.4-12 and 6.4-13 
Show the results obtained for uncontrolled and controlled responses for 
a positive and negative feed flow step respectively. In both cases the 
feedforward action can be seen in changes in steam and reflux flow 
rates. As can be seen the average offsets in bottoms composition are 
greatly decreased by the presence of the control action. The sustained 
oscillations evident in the simulation and experimental feedback tests 
are again evident. The effect of including feedforward control action 
should be to reduce the average offset in each composition and it can 
be seen by comparing these responses with those in Figures 6.4-9 and 
6.4-10 respectively that this is indeed the case, although the reduc- 


tions are small. 


6.5 Comparison With Previous Control Studies 

Berry [6] has carried out an experimental evaluation of two 
control systems on the column used in this study. One system involved 
the design of decoupler elements using a transfer function model 
derived from pulse data. The other system was a ratio controller which 
tended to remove the interactive effects by physical means. This ratio 


controller was a two point controller in that steam rate was used to 
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control bottoms composition and the ratio of reflux rate to overhead 
vapour rate was used to control overhead composition. 

Figures 6.5-1 and 6.5-2 (from thesis by Berry) compare these 
control systems for similar feed flow step disturbances. The feed 
Flow step magnitude was 15%. The controlled responses indicate similar 
behavior of these systems and both appear to be superior to the multi- 
variable control system investigated in this study, since both end 
compositions appear to be well controlled. This can be seen by com- 
paring the results in these figures with those presented in Figures 
O24 False. Los 

The reason for the superiority of the noninteracting control 
system is the inclusion of time delays in its design. The experimental 
open loop transients presented in Chapter IV indicate that there is a 
time delay between the initiation of a disturbance and the beginning of 
the composition response. These time delay terms were included in 
Berry's transfer function nicdel and appeared in the decouplers derived 
from this model. The original 20 state model exhibited these delays to 
some extent, since they appeared because of the series of first order 
terms. However, the reduced two state model does not have the number 
of terms required to introduce the delay to the same extent. This means 
that the control system designed from the two state model tended to 
over-compensate and yield oscillations, especially in the bottoms com- 
position control loop where the delay is largest. 

The ratio control system also gave better results than the 
multivariable control system. This is due to the inherent physical 
decoupling of this system. It tends to remove the effect of steam flow 


on overhead composition, which tends to upset the overhead control loop. 
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This allowed the loops to function almost independently as was the case 
for the noninteracting control system. The multivariable control system 
Should also account for the interaction but with the interaction time 


delays absent, the problem of cver-compensation and oscillations 


results. 
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CHAPTER VII 
DISCUSSION AND CONCLUSIONS 


The objective of this study was to obtain a process model based 
on transient balance equations and to apply this model to the design of 

a multivariable two-point control system. The 20 state model derived in 

Chapter III and tested in Chapter IV satisfies the requirement for a 

valid process model. It was found necessary to reduce the order of the 

model as outlined in Chapter V. The control strategy as outlined in 

Chapter VI was formulated using the reduced model and was used in sub- 

sequent experimental evaluation. 

The progression of study allows the following sequence of con- 
clusions to be drawn: 

1) The distillation column is a linear system with respect to product 
compositions only over a smal] operating region, that is for smal] 
input disturbances, so linearized perturbation techniques can be used 
in this region. 

2) It is not necessary to have detailed information regarding tray 
efficiencies as in perturbation form the slopes of the equilibrium 
data curves represent good approximation to tne slopes of the actual 
Operating curves. 

3) In the binary system used (methanol/water) the internal liquid flow 
dynamics cannot be neglected, and can be adequately represented by a 
first order expression. 

4) The distillation column used in this study can be adequately repre- 
sented by a twenty state model, with stage enthalpies and liquid 


flows as the states. 
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The two model reduction techniques tried in this study both yield 
acceptable reduced models. However, the two state model obtained 
by using Anderson's method was more attractive for the design of a 
control system as used in this study. 

The bottom composition measurement exhibited an appreciable amount 
of noise, due to the small composition of methanol and the approxi- 
mate peak area calculation. 

The combination of top product enthalpy (composition) control by 
reflux flow manipulation, and bottom product enthalpy (composition) 
control by steam flows manipulation yields a two point control 
system which exhibits interactive effects. 

The gas chromatograph separation and analysis cycle represents a 
measurement time delay. The application of an exponential filter 
to the noisy bottoms composition measurement has the effect of intro- 
ducing another time delay. The combination of these time delays 
causes a deterioration in the controlled behavior. 

It was necessary to modify the feedback control laws due to the 
bottoms composition noise. This feedback control system yields an 
oscillatory controlled response with a small offset in each composi- 
tion. This behavior has been predicted by simulation tests and has 
been subsequently verified experimentally. 

The inclusion of feedforward compensation with the feedback control 
system decreases the average offsets in compositions. The simula- 
tion and experimental tests verify this result, and indicate that 
sustained oscillations are present, as in the case of feedback 


control only. 


A comparison of this multivariable control scheme with those 
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investigated by Berry [6] indicates that the latter schemes are superior. 
This would imply that the method of model derivation, testing and reduc- 
tion and control law calculation yields a control] system which does not 
warrant the effort necessary. This multivariable approach is, however, 
the logical sequel to the simpler control schemes. This suggests that 
some method of process identification should be used to eliminate the 
initial steps of model derivation and verification. 

The major problems encountered in this study were time delays 
and measurement noise, both of which are due to the use of the gas 
chromatograph to analyze the bottoms. However, the properties of the 
dilute solution product are not conducive to other on-line analysis 


techniques such aS capacitance probe or refractometer. 


ton esob Afi matexe i 
eT9ySWON zi Toso WgQdsS ae eri aaa JAotts sit - _ 6M 


s7otg an .TSVSWOR .2uolsod ony scvlens OF Hoste 


me S03 norteluat es ia en bs nots 


2) erat ane ae yalont2 srt + isupge fetal sag 
wo! s 


7 


A ee ot 4, — 
fo ot healt eq bfucde Watsaitisnabi exs0074 Fo bortt am | Smo 


Ja2he19V Bak Notisvinab Tsbom to 2qak Tel 
. : 4 _ : 


i XbWte 2Tad nt borSsniegas. sel dove yotem silt 


to 920 SH7 oF SUD SB Hothw to Adocd ,seton fromeyyessm@ b 


f-no taddo ‘od ayfoubme ton sh touborq ootdulde 


TSTIMOTILATST 19 SGOVG S5Ns3 96GB. 2B Nave. Zann et 
= — 7 


: - 
= 


es 


CHAPTER VIII 
RECOMMENDATIONS 


As indicated in the previous chapter the process of model 
derivation and testing appears to be superfluous when the model is then 
reduced for control application. A more reasonable path is to derive 
a lower order model by the use of experimental data. Pulse testing has 
been applied to the column used in this study and a good transfer func- 
tion model obtained [6]. It is also expected that in the near future, 
pseudo-random binary sequence testing will be applied, and the model 
obtained will be compared to that derived from the pulse test data. 
However, if a model of the state space form is desired, it would pro- 
bably be much simpler to apply some form of learning or model-reference 
adaptive identification procedure. However, all these methods give only 
approximate models, when applied to a non-linear system such as a dis- 
tillation column, and the one chosen must depend on the operating range 
and application. 

The overhead product control configuration used in this study 
was to control the reflux flow and use the overhead product flow to con- 
tro] the condenser level. Recently some indication has been given that 
the reyerse scheme, i.e. controlling overhead product flow and using 
reflux flow to control condenser lJeyel, yields a more stable control 
system. It also has the effect of partially accounting for changes in 
overhead yapour flow, and so may produce a physical decoupling. This 
configuration should be tested and compared to the scheme used in this 
study to see if this improvement actually appears. 

The major problem encountered in this study was the noise and 


time delay in the bottom composition measurement. It may prove 
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advantageous to employ a Kalman filter or an observer to counter the 
effect of the noise. It may also be advantageous to use a predictor 
such as a Smith predictor which accounts for the time delays in the 
control system. However, both of these alternatives require the use 
of a good process model and so the model would need to be carefully 
evaluated before being applied. 

The non-linearity of the distillation system, coupled with 
the problems of process modelling, suggest the application of a 
learning or adaptive control system. This approach removes the need 
for a process model, and any control law derivation. Such a method 
could be applied providing some form of filter were used on the 
bottoms composition. 

Each of the alternatives to modelling and control presented 
above has been suggested as a possible improvement to the method used 
in this study. It is felt that each of these warrants some considera- 


tion and possible subsequent experimental study. 
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APPENDIX A 
NOMENCLATURE 


Chapter III: Model Derivation 


Scalars 
by sb. 
D 
F 
h 
H 
K 
L 


Xx 


y 


Vectors 


ies} 


te 


irs 


liquid/vapour enthalpy linearization factors 
distillate flow rate (1b/min) 
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liquid dynamics flow variation constant 
liquid flow rate from a stage (1b/min) 
heat loss from a stage (BTU/min) 
condenser heat duty (BTU/min) 

reboiler heat duty (BTU/min) 

time 

temperature 

liquid flow dynamics time constant 
vapour flow rate from a stage (1b/min) 
stage mass holdup 

liquid phase composition (wt % MeOH) 
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Subscripts 
1-10 
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G 
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IE 


Superscripts 
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Chapter IV: Model Evaluation 


Scalars 
a - scalar example coefficient 
fo) - constant 
b - scalar example coefficient 
B = constant 
d - scalar example input disturbance 
g - scalar example state coefficient 
h - scalar example input coefficient 
k - number of incremental time intervals in 
model time base 
r - eigenvalue of matrix G 
= tlme 
i - model time base 
Xx - scalar example state 
Vectors 
Xx =6.Gabe; VeCto), 
D - disturbance/control vector 
z - canonical state yector 
Matrices 
G - continuous model state transition matrix 
H - continuous model input weighting matrix 
M - similarity transformation matrix 
‘ - discrete model state transition matrix 
: - Jordan canonical form input weighting matrix 
$ - discrete model input weighting matrix 
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Chapter V: Model Reduction 
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exponential filter constant 


designation of general time interval 
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vector of control variables 
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continuous model disturbance weighting matrix 
discrete model disturbance weighting matrix 
feedback control matrix 

feedforward control matrix 

setpoint control matrix 


discrete model state transition matrix 


Laer 


“notsey Sietebsvariti~ «gk 
(orev sdeye btu 265N : i on 
2 


ey ! 
ayrfroarse atade Po tar9ev = 


xivtam fortteasys ‘Sopa2 lgbometouniines - 


xritem potidetew 2forns Fabom -euournisi09 
Arian ehisipiaw 2faysnos (spon atsaaio : 
xiysem pnitdptaw soneqyt eh ishow Swouninus: ~ 
aMdpm pnitngisw sansdiwserb: fabom sdsypett - 
arvtam Jowtwes A06dbsaT - 

Kivjam lorde> biowtotbest - 

afitem Toritioa sriogise - 

xiryiem nots fanert jagsce ishom.staroetb - 


ag Oe ee “24 


135 


APPENDIX B 
MATRIX REPRESENTATION OF BALANCE EQUATIONS 


The dynamic equations outlined in Chapter III represent two 
sets of ten coupled differential equations. These can be written as 


two matrix-vector equations as outlined in the following sections. 


B-1 Overall Balances 

The total mass balance relations represent ten differential 
equations relating the changes in mass holdup in each stage to the 
Stream flows to and from that stage. The following is a summary of 
these equations: 


a) Reboiler (Stage 1) 
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c) Feed Tray (Stage 5) 

Wes + Wie = Va + Le - Vp (es 
d) Condenser (Stage 10) 

Weide Mies Gu yoo 


These equations can be written as a matrix-vector equation of 


order 10 of tne form: 


AB 70 MONTATHReSAGaA x 7" 7 
a ye rT 7 _ 7 
p sat 


awd Treas Be ett. 2nodivewps oF T 


_ 


26 i993 Fv 9d 62 sdaut «enobteups Tatinovatt th .bs aac: ES Ye efee 


.andidss2 patwor to? any nt Bontistuo 26 20} $sup NORDSY-Hh sm OW: 


ssonstid a i 

[stinsvsttth nad snoeews4 enobdelay Sonetad czam Tedod sit 7 _ 
ot ot spate dose nt qublod ezem nf 2eprcilo av paitaatan ar > Jaupa 
to yyseme.5 ef ontwolfot aff “spss2 Jentd wort BAS GF aio ati 

‘ 2noTdsups 4 - od 

({ gps72) st fodah f 


i ! : re i 7 = 


. 7 Y & = } i 
, = » ° al al t fo" 


(ph apie) yey! leq 91 y Aa 


aah yr + VS swt on 
Po y Cee riz! rj" + Faw 


; ae 
(2 spate) yerT b ail te 


: : 
‘ i ' ' ie 
Pt gla he att = a + oe 
: uf 7 
‘ ; : Of spsd2) vaansbne 
i n 7 7 ¢ — m 


+ 
yi iP 1 i 
fm pi = oV os 3 vo 
ss ial act ei ad ne se ae sea i 


a yi « ot eee i ot y 
' : a) Sr »” 7 = . 


=m -N -M--T -N -00 - BR -O - Oo - 
(E31 cy, een ener ae eg ia 


0 
0 
0 
0 
1 
2 


er icf te) ©) ©) ! 
SS BOO RO Og IO eS 
OF Ow, OPO ae Te BO 
Sa ZF oS a ea = a (a Vinay = er a =) 
SS) SS) Be Se Cc ea © 
Cy OS = =O $O je is) © 
a & = Es es Se &| @& @& & 
ME hess ey Te) j~) (Ss) Clue) eS) 
errr ate @2 & @& © © © 

(ey te) (ey ee) eer me fe tea) = 


ey Ve) Key Ke) fe) Sy Key (se) @P XS) 
JF Se oe oOo 2 S| 2 & &S © 
J Oe 2 E& S&S © OS S&S ©€ 
J er FQ & & @ © So ce 
SehOU SOT R=) Sane eee see 
se eo —= ©& © & S&S © 


Ste SIN) Sinn) Sear Spl SWe) SS See) Sepys 
S SoS SS SSS SSS SS 


[atk Cee Sk Pe eo ee EE a eG 


losiaio Rouge non ohay Al 
se © ©& © © © © es 
Oe eS eS 
OF fe OE ON IO ar 
SE ee eS ee ag 
CES (ee = 94> 1S 
Bg OS ES OS 
OP Oe Ss OG Oe > 
OF ee ee eee 
Ss = ©) © oo) @) OF (© 


136 


TOTAL iWASS BALANCE 


Figure B=. 
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This equation is shown in its expanded form in Figure B-1.1 where A, B 
and C, are defined. The vector D includes terms not appearing in the 
total mass balances but which appear in the enthalpy balances and are 


either controls or input disturbances. 


B-2 Enthalpy Balances 

The ten enthalpy balance equations represent a set of ten 
equations, similar to those for the total mass balances. They give the 
rate of change of enthalpy content in each stage as a function of the 


stream enthalpies. The following is a summary of these equations: 


a) Reboiler (Stage 1) 


b) General Tray (Stage I) 
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d) Condenser (Stage 10) 
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These equations can be expressed as a matrix-vector equation of order 10. 


The equation is as follows: 
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The matrices and vectors in this equation are defined in Figure (B-2.1) 
and the nomenclature appears in Appendix A. The term hio which appears 
in the vector D, allows equation (B-2.1) to be written in the same form 


as equation (B-1.1). 
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FIGURE B-2.1: 


ENTHALPY BALANCE MATRICES AND VECTORS 
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C-] Experimental Control Configuration 

The distillation column used in this study is completely 
instrumented with local analog controllers and data display recorders. 
It is also interfaced to an IBM 1800 Data Acquisition and Control com- 
puter which can be used in an analog supervisory mode for control and 
data gathering purposes. A schematic diagram of the distillation 
column showing the complete instrumentation is presented in Figures 
A, B and C with an accompanying legend. From these figures it can be 
seen that numerous different control configurations are possible. The 
local analog controllers were all operated in the proportional-integral 
mode. The controller settings used are presented in Table C-1.1, and 


were found by on-line controller tuning. 


C-2 Sample Experimental Data 
C-2.1 Steady State Data 

As indicated previously, the distillation column used in this 
study is interfaced to an IBM 1800 Data Acquisition and Control com- 
puter, to enable continuous data gathering. A computer program has 
been written by previous workers to sample all measurable variables at 
a given interval and do a steady state analysis of the data. A mass 
and energy balance is performed and closure errors and heat losses cal- 
culated. A sample of the output from this program is given in Table 
C-2) |. 
C-2.2 Steady State Parameters 

In addition to the data collected, several additional para- 
meters of the column operation are necessary for the model. Utilizing 
all the data collected coupled with a knowledge of the physical pro- 


perties of the system, it is possible to calculate the values of 


ma res 
“m2 fortnod. bn notte we fr oo 


org fortnoo 18k, abom ryOe 


a Td 


afl -.sfdi2209 31s eioiSe eee rs snoichie hone 
fsreotat<lenatiroqorg sont mf bate raqo (16 arew 2ya{ Joven: pOTenE teat <5 
bone ./ c1-0 sfdel. nh tedno2er9q 6X6 bSa2u conesbe! yal i et7ng, SAT ‘ebom 


-pitrd set iditrioa anil=ne- i ae 


“east ner aS 
sethpigte ves fo 7 


ba HU 


“4 


. 


- , 


Silt ni? be2u omuloo nottel hfierh sit vi eporvedg bad 


$6 @hlcstisy ol Usrvvesom [la al amibe OF 2yaNvew . eng t¥ 


Zech A> .oisb Sit fn zreyiam state ybbat2 6° ob bye Ts 
eas fae 7 = 


+(go 2s@20! J8$h bre ey074s sviedts bas Bonnoived 24 tee 
Hfapl ni nsvip 2iompwpesy eid jrovt dugsao sat to orang Le 


Sad 


me 
* 


Miia 


sara: Hirai Fsvavae ,betost {oo steb sien 
ae sem nies own or 


141 


NOILDIS GVAHYSAO - NWN1O0D 
| NOILVT1HLSIG JO WVAOVIG DILVW3HDS ‘V 3FYNOlS 


© CeOne, 
8 
6u:6607 & 
4@SUSPUOD JUBA 
; GUO Jara] 0}) @é2u4 i 
AN3A jO4}U0D yUIO 


H 38S 10 6ulHGo4 (2 3 
a ay © 


SASS ORES SSne 


AS? i} 
ong 
H ulesg 
Ole | 2 
Me speteonpee cenacaresetesests | 
| Fears © 
(i) Gi) Za “ wy 


eeuel. 4. bi ok Pi hag Ronioes? 
|_| aaa CY) (tz ~ =e a wan | ot e}e) 
Ly “33Y o~ REY, Q 6HulhGo04q 


yu 


> 


\* a7 


‘ 
. 


a 


rer 


- 


4 


fea to gniggol 
— Jortnes tniog 
ong level of) $3H7 


i 
i 
: 
. 
‘ 
4 


(reenebroo tnev : 
' pniggos 
: x SAT —~ 
oF <2 


ieee) 
er 
ae as ee ee ee ee Oe 


VOITAIIUT2ZIO IO MAROAIG DITAMSHIDe A FAUOIS 
UOITDI32 GAZHAIVO - AMUIOD 


142 


NOILD39S -dIW - NWN105D 
NOIVTINILSIG JO WVYOVIG DILVWSHDS 9 ‘€@ FINO 


(auo jaAa] SHUR} 0}) or eR aera 
1onaoud dol /ha\ wT) hh 
{1 ess 
Ses i ! 
O--tw 2 LXH Die 
, Las 
| -G q334 mz 
| 8 ERY 
| : 
i zOu4 
| mitt — os 9 
Te, Pes 
, ob | | Oe xX 
, | AH aac eb 
| to Xx xX (iene 
: (4) Oe 7 uleig (9) ey ® 
uleig = 
| aes 4 <O 
eis [pe "75a 70]= Bar a (x) | ey 
| £XH ae i ae 
| Sie a 1 a eo 
| 


lace Ange ateems W nar q ad ' gear! 
: i ; 4 eB he 
TaAvOeR= F075 a Pet, Oe wed 
| afar eye) CArent at y ae es . Wise 


VOWIRINT2°G 9S MAROAKI SITAMGHD? -.6 Rug 


ADT 332 - FIM - AMUIGS 


143 


NOM aSewWOle@d= NWiloS 
NOUV IMISIO SOuwWvdoOVIGSSUVWaHOsSS Sado 
(b Janay \ pat) 


Bulbbo; Sf 
uolzepnojes 
Hulbbo; 
|| S4ue} 0}) 
19ndqo0ud ' (i4) 
WOLLOG 
Se 


uleig 
uleig ee Hulbbo; 
yalVM 
; OINITOO9 
\ 
= s ajesuapuoos NO || eee Le ee abund 


Ged eens (ee sosma 


\ ay 


\ \ 
ee en On | ae oa a 


v : 
‘ 
tea 
a6ind | “ y [aAa], 494 
r tee 
4o HulHhbo| Old) 


E 
Es 


. VOITALT2T IO MAR AIG DITAMBHOS: <3 SIONS 
MOTD 72). MOTTIOS > HMBFIOD 


1h 


YIINdWOD OL SANIT NOISSIWSNYYL YiVO =—=———— 


IN ERORES BSc OQ CA<__A< ee SINI1 SS390Ud YILVM ONII00D =o oc ome 


JONTBHAG IOS ICI. SIN] 1 LNIWNYILSNI S$S3900Ud 


SAINIMES S29 Oda dad My aON Via Wi ——————oeneeeee SSINIU AT LU ENAUIARSIIMS ING OMA SYM) Se 


SINI1 $S390U8d 


YALNdWOD AB (U-UU FIGVIYVA) TOYLNOD GNY 9NI9D901 dWNd 


YILNdWOD Ag (U-UU FIGVIYVA) ONIDIO YJINGSNVYL 13A31 


JAIWA QIONITOS YITIOYLNOD/YOLVIOIGNI 13A31 


Ch) e) & 


YITIOYINOD/YICUYOIIY JANLVeIdWw4l 


> 
a, 


JAIVA GNVH 


YIddOIIY JaNLVAIdWILl INIOdILINW (1182 d/P) YJONGSNWYL MOIS 


M01 


YIINAGSNVAL 34NSSIUd YITIOYLNOD/ssdyODIY MO14 


is 
(oa) 
36 


YIINASNVYL JYNSSIYd Y¥id¥yO9I4¥ MOT4 


2060808 


Ded & €) eb © 


HILIMS FYNSSAYd Sd (aqgoud aduejiseded) YIONGSNVYL NOILISOdWOD 
JAIVA ONILVINDIY JYNSSINd Add (a9duey9eded) YITIONINOI/YIGYOIIY NOILISOdWO9d 
YFTIOYLNOD/YOLVIIGN! JYNSS3IYd id (ydesboyewosyd) YIGHYOIIY NOILISOdWOd 


-GN99391 


~*~ _ 


A208 TAOIROTADIONY TAUZ2IR4 & : 


rmerganensd rs ng 
_ 


= GVIAV DRLTAIUOAY WUSTINGA VHA teanehaeqgesa! RILIGRIVOD 910 0238-401T 12 ORO: cy 


— ; GGTttWe aeuzeteq = ¢8 iadote 427Etl 2eg2 2? VITUTP YAS voineaawne 78) 7 


he. 
oe a : 2 
VIDOE TevARE [Hes TAG 4 : aieRooad MOST + SS - 
a / -_ 
Hom at aor ks rah eq liam ) S340 mos RAOROIIE work | ey 
wa ‘ ae ~ : _ , 
iba 7 - 
Ga eR ABUTARSA MATUTHI OSE IUM a 4 /iad Gib) SIZUGEVART YL039 e, at 
j ee — 
z al : ee . | 


= | RIGOR HO DRIER DIF TAL IAWAI MIT 1e 4) WIAVRV RB “WES 


: inne \ ~ 
7 ; - WJAY GiGekgj0z We } way O@MMaIMOTADIONT Javad 


: 
fe 


- : ms bs 
7 AaTUeMOD Y= in-nm AISATHAVA OMI gO0s FAS RFD UGOMAAT BY D4 (B) 


7 
- 


i =n 8 
a AIUTMUD Fa (n-An JIEAI FS JORTMOD CHACDMI deGd frre : ‘MOS > A-49 > 


- 
a 
e2720KS 
vas 
Tit. 
a 23 ; RITAW iAP i" — a Csyvil ‘ Ay As 11 TAMUSAS — a 
vw! €<3 njirt 
7 
: | * _— 
" i q — 
V4 on eee ee ee a3 ish evry ,D0R4 — 
4 
5 2 » Oo2T See irr 
i Lan 34 = c. a Wij 2 Lé a rAW SHIJUUL oo + oe soe " 
R27 ow aa) eat { PUM, ART TAG —_ — => =o 
i 
g = 
— 


TABLE C-1.1: 


Controller 

Feed Flow 
Reflux Flow 
Steam Flow 
Column Pressure 
Reboiler Level 
Condenser Level 
Feed Temp. 
Reflux Temp. 


Overhead Comp. 


TYPICAL CONTROLLER SETTINGS 


Prop. Action 
Kp 


Wie) 
100 
Zo 
50 
50 
50 
oo ee 
Pet) 
-200 
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Int. Action 
(reset time) 
(min) 
0.8 
0.4 
0.1 
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variables which cannot be measured such as internal column flowrates 
and compositions. However, difficulties can arise in estimation of 
tray efficiencies and equilibrium factors. In this study the stage 
liquids were sampled and the samples analyzed by use of a gas chroma- 
tograph. This composition data, with the assumption of saturated phase 
enthalpy relations, allowed simple calculation of values for all vari- 
ables. 

The steady state material and energy balance equations for 


a general stage are: 


Ly Mae = Ee Nats. (C-1) 
Ba i Seay Vip set a 
EO = aa ae ey ies (C-3) 


where S. is a sidestream input to stage i and Qi; is the heat loss from 
Stage i. The heat losses were obtained from an overall neat balance 
and were distributed as indicated by Svrcek [46]. 

With the measured liquid compositions and heat losses, and the 
enthalpy relations, there are five unknowns in equations (C-1,2,3). 


on stage) 


However, if values of Ley and Vv. (from the stage above the i 
were available, the equations could be solved. Since a steady state 
balance on the condenser yields values for Vg and Lige the equations 
that apply to the top tray can be solved, and subsequently each tray 

can be solved, working from the top of the column. A computer program 
was written to employ the available data to calculate the internal flows 


and compositions. A sample of the output from this program using the 


steady state data given in Table C-2.1 is presented in Table C-2.2. 
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TABLE C=2.1: 


SAMPLE STEADY STATE DATA 


STEADY STATE DATA 
RUN NO SAMPLE 
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TOSS 72 
FEED FLOW 2e135 LB/MIN BOTTOM PROD Pecb2  EBMIN 
REFLUX FLOW 16956 LB/MIN TOP PROD 1.002 LB/MIN 
STEAM FLOW 12827 LB/MIN COOL WATER 68.988 LB/MIN 
Pema PAT e 4 FEED COME 45249. WTOP C 
TOP PROD 96.88 WT P C BOTTOMS COMP Les Op INT wPeaG 
BEED UNC ET 163.4 DEG F REFLUX INLET Paired, DEG SE 
STEAM TEMP 220.69. DEG. EF PRESSURE O.7 IN H20 

Mick. Tweet OAL BUAUL, A IND CE 
FLOW COMP METHANOL WATER 

(LB/MIN) (WT PCT) (LB/MIN) (LB/MIN) 
Riset) Zeiso 452499 0.971 Peo 
BOTTOM PRODUCT 1.262 1.350 0.017 1.245 
LOPLPRODUCT 1.002 962889 ORR 0.031 
CLOSURE ERROR=PC 6.0 1.6 Oreall 


EIN Ek Guay B AL, ALN CE 


ENTHALPY IN ENTHALPY OUT 


(BTU/MIN) (BTU/MIN) 
COOLING WATER 3733.3 5058.7 
REFLUX 189.8 203-2 
TOP PRODUCT 104.1 
FEED 308.3 
STEAM 2184.6 418.2 
BOTTOM PRODUCT 264.3 
TOTAL 6416.1 6048.7 
HEATSEOSsS 367.4 
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TABLE C-2.1 Continued 
STEADY STATE CONDITIONS BASED ON 20 POINTS 
RUN NO SAMPLE 30/3/72 

FEED FLOW = 26231 \LB/MIN DEV= 0.0041 
REFLUX FLOW = 2.233 LB/MIN DEV= 0.0093 
STEAM FLOW = 1.2827 LB/MIN DEV= 0.0082 
BOTTOM PROD = 1.264 LB/MIN DEV= 0.0823 
TOP PROD = 1.144 LB/MIN DEV= OS0L17 
COOL WATER = 68.919 LB/MIN DEV= 0.0641 
TOP PROD = 96.889 WIPC DEV= 0.0811 
BOTTOMS COMP= 1.350 wWrPC DEV= 0.0000 
FEED COMP = 45.499 WIPC DEV= 0.0000 
PRESSURE = 0.758 IN H20 DEV= 0.0968 
COND LEVEL = 5.055 PSIG DEV= 0.0246 
REB'R LEVEL = 9.759 PSIG DEV=" 0.1180 
DIFF PRESS = 10.458 PSIG DEV= 0.1520 
REBOILER TEM= 208.3 DEG F DEVS \Oes 700 
PLATE 1 TEMP= 193.1 DEG F DEV =O rg. 
PLATE 2 TEMP= 178.9 DEG F DEV= 0.3384 
PLATE 3 TEMP= 163.6 DEG F DEV= 0.2786 
PLATE 4 TEMP= 166.63 DEG F DEV= 0.3099 
PLATE 5 TEMP= 158.0 DEG F DEV= 0.2833 
PLATE 6 TEMP= 152.8 DEG F DEV= 0.2603 
PLATE 7 TEMP= 150.0 DEG F DEV= 0.2653 
PLATE 8 TEMP= 147.5 DEG F DEV= 0.2690 
COND TEMP = 144.5 DEG F DEV= 0.3314 
STEAM TEMP = 228.9 DEG F DEV= 0.4987 
COND'T TEMP = 226.6 DEG F DEV = 2032 44 
REFLUX FLOW = 121.7 DEG F DEV= 0.2714 
FEED FLOW = 89.6 DEG F DEV= 0.2161 
BOTTOM FLOW = 99.0 DEG F DEV= 0.1767 
REB O'HEAD = 206.6 DEG F DEV= 0.3682 
FEED INLET = 163.4 DEG F DEV= 0.2947 
REFLUX INLET= 137.2 DEG F DEV= 0.4435 
COL O'HEAD = 148.2 DEG F DEV= 0.2935 
WATER INLET = 54.1 DEG F DEV= 0.2801 
WATER OUTLET= 73.4 DEG F DEV= 0.3059 
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TABLENC-2.)2tee@SAMPLE STEADY STATE PARAMETERS 


DISTILLATION COLUMN 
STEADYSTATE DATA 


FLOW COMP 

(LB/ MIN) WT FR) 

FEED rae Me) 0.4550 

DIS VTECATE 1.002 0.9662 

BOTTOMS Pes 0.0135 

STAGE LIQUID VAPOUR 
FLOW COMP FLOW C OMP 

1 Lel337 0680135 Zel52 “O.0579 
2 36285 (020390 ZecOt "O01 52 
3 32340 0.0970 22-346 0.3619 
4 36479 022450 20429 044957 
5 32562 0.3390 Zee 20669715 
6 1.650 0.5345 2-838 0.8052 
if, 16836 OO. 7175 20360 0.8633 
8 1.958 0.8410 2-984 0.9321 
9 1.982 0.9150 20978 0.9662 
10 1.956 0.9662 0.000 0.0000 
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C-3 Sample Model Parameters 
C-3.1 Liquid and Vapour Volumes 

The model as derived earlier has as parameters the steady 
State liquid mass holdups. Values were calculated from a knowledge 
of the liquid volume in each stage employing a composition dependent 
density relation. The condenser and reboiler volumes were found by 
filling each with liquid to the controlled setpoint and measuring 
the volume of liquid by draining each vessel. The tray volumes were 
more difficult to obtain as observation indicated that there was more 
liquid (and foam) on the trays adjacent to, and including the feed 
tray than on tne others. A weighting factor for each tray was esti- 
mated by observing the holdups and was used to compensate for varia- 
tions in liquid loading throughout the column. Tne procedure followed 
in determining the tray holdups was to operate the column at steady 
state and then abruptly (and simultaneously) stop all inputs and out- 
puts. The column was then drained and the total tray volume deter- 
mined by subtracting the previously measured reboiler and condenser 
volumes. The individual tray volumes were then arrived at by applying 
tne weighting factors. Table C-3.1 presents the weighting factors used 
and the volumes allocated to each stage. 

The vapour yolumes were also required, to estimate the vapour 
enthalpy contribution to the transient heat balance. These volumes were 
found by calculating the total volume between stages from a knowledge of 
the column dimensions and then subtracting the liquid volume holdup. 
Table C-3.] also contains the calculated vapour volumes of each stage. 
C-3.2 Enthalpy Linearization Factors 


As outlined previously, in Section 3.5, the vapour enthalpy 
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TABLE C-3.1: LIQUID AND VAPOUR STAGE VOLUMES 


Liquid Volume Liquid Volume Vapour Volume 
Stage Weighting Factor (ft3) (ft3) 
] - .2474 Jagec 
2 1200 .0259 ©4151 
3 1-05 .027] .4139 
4 Latte) .0285 24125 
5 1S .0298 4112 
6 120 0311 .4099 
7 1.20 “0311 4099 
8 Pos .0298 -4112 
9 1210 .0285 1.090 


10 = 20557 12302 
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perturbations can be expressed in terms of the liquid perturbations by 


an expression of the form: 


Ho = by h 


where by is the slope of some pseudo-equilibrium/enthalpy relation for 
each stage at the initial steady state values. 


In general 
Te eta) 


hiereee “f X%s) 


2 


and more specifically for constant pressure and stage efficiency 


1 nd 

xX = Fy (T) 
so that: 

eae ae iety) 

h. = f, (x) 


So the factor by can be expressed as: 


b, = dH = ES iia es 
Neos, dy dx dx S68 


The physical property data for the methanol-water system have 
been presented by Svrcek. Although the phases in each stage are not 
in thermodynamic equilibrium, if a constant efficiency is assumed, the 


slope of the actual operating curves will be the same as the equilibrium 
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curves. So the equilibrium data can be used, and the assumption of con- 
stant efficiency should introduce only a small error. 
Suppose that stage 9 is considered. Table C-2.2 gives the 


compositions to be: 


96.62 wt% 


< 
It 


91.50 wt% 


x 


so from the data presented by Svrcek; 


a = - 613.8 (linear over entire range) 
a = - 10.0 (linear between nearest tabulated values) 
we = .342 (linear between nearest tabulated values) 
therefore 
] dy dxa" dx 


_ (-613.8)(.342)  _ 
<= aa 20.95 


This yalue and the values for the other stages are given in Table C-3.2. 
All values were calculated using the compositions presented in Table 
C-2.2. It can be seen from Table C-3.2 that there is a large variation 
in the parameter by from stage to stage. This is due to the S-shaped 
y-xX and h-x curves. 

C-3.3 Liquid Dynamic Constants 


The model outlined in Chapter III of this work accounts for 
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TABLE C-3.2: ENTHALPY LINEARIZATION FACTORS 


Stage Dy 
] Lae Peso, 
2 19.40 
6 Sesit) 
4 Ne We 
5 2405 
6 8.94 
7 8.38 
8 18.62 
9 20295 
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liquid dynamics by the use of first order expressions of the form: 


where T. is a flow change time constant and K is a final steady state 
flow correction factor. The K; term is necessary since, in general, 
the liquid flow from a stage will not change by the same amount as the 
flow onto the stage for input changes. The estimation of these para- 
meters t. and K requires data for the internal liquid flows as a func- 
tion of time. This type of data, for the column used in this study has 
been given by Svreek , and so was used in the estimation of the 
required values. Table C-3.3 presents the calculated values which are 
averages of values determined for several sets of data. The large value 
at the feed stage is due to slight subcooling of the feed. 

The feed tray liquid dynamic equation has one extra term for 


the feed input and is of the form: 


et a ike i} { * L 
bs = ay Ete Ks * F Kes a 
Also the distillate flow is dependent on the feed and reflux flows and 
the expression is: 

et ] | I i} ( 
De ee i” Sipeep Era 
The yalues of the parameters in these expressions were also estimated 


from Svrcek's data and are also given in Table C-3.3. 
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TABLE C=3.3: 


LIQUID DYNAMIC CONSTANTS 
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C-4 Reproducibility 

Some investigation was done to establish the reproducibility 
of steady state conditions, and of transient data. The steady state 
reproducibility entailed the definition of an initial steady state and 
Subsequent introduction of a pulse disturbance of about 10 minutes 
duration. The column was then allowed to come to steady state. The 
overhead composition was used to indicate the steady state reproduc- 
ibility. Table C-4.1 presents a set of data taken for four consecutive 
disturbances and the type and duration of pulses used. 

From tnis table, it can be seen that the deviation in overhead com- 
position can be as large as .3 wt % MeOH, but a deviation of about .] 
wt % MeOH would be average. 

The reproducibility of a transient run was investigated by 
comparing similar open loop cases. Figure C-4.1 compares two similar 
transient runs. As can be seen the overhead composttions at final 
steady state are different by about .15% and the bottoms composition 
differ by 0.5%. From these observations and consideration of the data 
in Table C-4.] it can be concluded that the transient data is about as 
reproducible as the steady state data, having variations of about the 
same order of magnitude. 

The reproducibility of a multivariable feedback control run 
was also investigated. Figure C-4.2 compares two such runs. Although 
the initial steady states are different, the controlled behaviour is 
very similar. The bottoms composition shows similar oscillations with 
an average offset of about 0.5%. The overhead composition in each case 
has very little offset and small fluctuations. This would indicate 


that the controlled behaviour is reproducible despite slight variations 
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in steady state. 

The main difficulty in attempting to reproduce experimental 
data is in the slight variations in conditions which lead to differing 
results. Slight differences in tray loading or in steam pressure can 
Cause appreciable differences in transient or final steady state read- 
ings. These problems however do not appear to introduce appreciable 
problems into the control studies, as controlled behaviour is fairly 


reproducible. 


C-5 Model Matrices 
The twenty state model, in differential form can be written 


as: 


The coefficient matrices for this formulation are presented in Table 
C-5.1, and were calculated from the steady state data previously 


presented. 


C-6 Noise for Simulation Studies 

The simulation studies presented in Chapter VI of this work 
include experimentally measured noise in the bottoms composition. This 
noise was determined at an experimental steady state and the case used 


is presented in Figure C-6-1. 
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TABLE C-5.1 MODEL COEFFICIENT MATRICES 


MODEL OX/DT = AXxX + BxU 
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APPENDIX D 
COMPUTER PROGRAMS 


The programs written and implemented on the IBM 1800 Data 
Acquisition and Control computer during the course of this study fall 
into five general categories: 


1) Standard DDC operation and associated utility programs. 


) 
2) Steady state data analysis programs. 
3) Test run data accumulation and control programs. 
4) Model derivation and reduction programs. 


5) Miscellaneous utility programs. 
These programs are briefly described in the following sections, and these 
sections have been combined with program listings as a Users' Manual. 
Table D-1 presents a tabulation of the programs and data files and their 


uses. 


D-1 DDC Utility Programs 

The system DDC (Direct Digital Control) program is interfaced 
to the distillation column so that many variables can be displayed and 
the feed, steam and reflux flows controlled. Table D-2 is a tabulation 
of the DDC loop designations. The data accumulation loops are turned 
on and off by the program RGM40. The data collected is punched to cards 
with the titles stored in RGMO1 by the program RGM41. The program RGM42 
is similar to RGM40, but is used to turn the data acquisition and 


control loops on or off. 


D-2 Steady State Programs 


In this category of programs there are two sets of programs, 
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DISTILLATION COLUMN PROGRAMS AND DATA FILES 


MAINLINE CORELOADS 


BALNC= CORELOAD TO DO MASS AND ENERGY BALANCES AND OUTPUT 


DASS = 


CVIEUED CORELOAD ‘FOR STEADY STATE (DATA *COLIZECT ION 


DATAC= CORELOAD TO GATHER STEADY STATE DATA AND WRITE IT ON 


DIS iis SPANO Dp. 


DCMNAD »OCMD1,DCMD2—- BATCH CORELOADS TO CALCULATE MODEL MATRICES 


DGRED= 


Chiknt h= 


GRASH= 


RGM40—- 


RGM41—- 


RGM42- 


RGM43- 


RGM45—- 


RGM46—- 


GM46A- 


RGM4A- 


RGM47- 


RGM50—- 


RGM51- 


BATCH CORELOAD TO PERFORM LEAST SQUARES REDUCTION 


INTERRUPT *CORELOAD FOR GeCe OPERATION 


SYSTEM PROGRAM FOR EMERGENCY SHUTDOWN 

TURN ON/OFF DDC DATA ACCUMULATION LOOPS 

PUNCH DATA FROM DDC DATA ACCUMULATION LOOPS TO CARDS 
TURN ON/OFF DDC DATA ACQUISITION AND CONTROL LOOPS 


TURN ON/OFF EMERGENCY SHUTDOWN ECO 


START AND STOP BOTTOMS GAS CHROMATOGRAPH 


TRANSFER DATA FROM DISK FILE *"RGMO3* TG *RGMOs* 


PUNCH DATA FROM #RGMO4* TO CARDS 


READ LABELS FOR GM46A INTO 'RGMOL! 
TURN ON BOTTOMS G.Ce FOR ONE ANALYSIS CYCLE 


INTERRUPT CORELODAD FOR CONTROL AND DATA ACCUMULATION 


BATCH CORELOAD TO READ CONTROL MATRICES INTO 'RGMO3B! 
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DISTILLATION COLUMN DDC LOOPS 


DATA ACQUISITION AND CONTROL LOOPS 


0601 
0602 
0603 
0604 
0605 
0606 
0607 
0610 
0611 
0612 
0613 
0614 
O65 
0616 
0617 
0618 
0619 
0620 
0621 
0622 
0623 


PEED: FLOW 

REFLUX FLOW 

STEAM FLOW 

OVERHEAD COMPOSITION 

TOP PRODUCT FEL OW 

BOTTOM PRODUCT FLOW 

BOTTOM PRODUCT COMPOSITION 
REBOILER TEMPERATURE 


TRAY LL TEMPERATURE 
TRAYe 2 TEMPERATURE 
TRAY 3° TEMPERATURE 
TRAY 4 TEMPERATURE 
TRAYS TEMPERATURE 
TRAY 6 TEMPERATURE 
TRAY Too TEMPERATURE 
TRAY 6 TEMPERATURE 


CONDENSER TEMPERATURE 

REFLUX TEMPERATURE 

FEED TEMPERATURE 

COOLING WATER INLET TEMPERATURE 
COOLING WATER OUTLET TEMPERATURE 


DATA ACCUMULATION LOOPS 


O20 
0202 
0210 
Ozeist 
O22 
0213 
O21 
O215 
0216 
O20, 
0218 
O2Z19 
0220 
O22) 
Geez 2 
0223 
O22Er 
0225 
0226 
0227 
0228 


GeCw DATA ACQUISITION 
GeCe DATA ACCUMULATION 
FEED FLOW 

REFLUX FLOW 

STEAM: PLOW 

TOP PRODUCT FLOW 
BOTTOM PRODUCT FLOW 
TOP PRODUCT COMPOSITION 
REBOILER TEMPERATURE 
BOTTOM PRODUCT COMPOSITION 
FEED TEMPERATURE 
REFLUX TEMPERATURE 
TRAY Cl TEMPERATURE 
TRAY oc TEMPERATURE 
TRAY 3 TEMPERATURE 
TRAY 4 TEMPERATURE 
PRAY, SatvEMPERAT URE 
TRAY 6 TEMPERATURE 
TRAY (7 (TEMPERATURE 
TRAY 8 TEMPERATURE 
CONDENSRR TEMPERATURE 
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one which reads the data from the column at steady state, and the other 
which uses the steady state data to calculate all unmeasurable variables. 
The first function is performed by the programs DASS, DATAC and BALNC. 
The first is a teletype queued coreload and is used to enter certain 
information and start the data gathering. Adequate description is 
printed on the teletype before each input request. DASS queues DATAC 
wnich gathers the required data for a predetermined length of time and 
writes the data to disk file PACOl. When enough data has been gathered 
DATAC queues BALNC which performs mass and energy balances using the 
data in the disk file and outputs the balances and variable averages. 
The programs are queued to V-Core and output in to the line printer. 
The second function outlined above is performed by the batch 
process program SSDAT. It does an iterative mass and energy balance 
on each stage given the input liquid flow, the output vapour flow and 
the liquid composition. This allows each stage to be solved in turn 
from the top of the column. The subroutine REGFL uses a false position 


search to solve the stage equations. 


D-3 Data Accumulation and Control Programs 

There are many programs in this category and each will be 
mentioned in turn, in approximately the order in which they would be 
used during a test run. The hardware timer which repetitively queues 
the accumulation and control program, RGM50, is started and stopped by 
RGH52. The type of control and steady state data are also set by this 
program. The control matrices are read into the disk file RGMO3 by 
the program RGM51. The model order, feedback, feedforward and integral 


control matrices are supplied to the program. The steady state data 


wr eee Ye 
7 7 7 - _ Hi 
7 i 7 
sorido sits: brs voted beste Jp amufoo oat mor? dab eat abet tate 
.2atdnitsy ofdewwesanm: Tis Stsfualss of ersb stad2 ybaste ett su 
.OWJAa bos DATA .22A0 2m 1BONG ‘ait vd bomriotieq 2t mofsonut tevit 2 
nteinoo 19tn8° 0d beew 2t bné beofso> boveup sqysefat 6 2f tert aap 


= 
2} notiqtiozeb steupsbA  .painsiitsp stab ony tute bas norsaarvotnl aT) 


oe 


JATAG 2ousup 22Ad .d25upey Suyn? Wasa sited sqysels? |a3 Ho bstnivqg - / 
bas omrd to dipusf banrerisiadsiq 5 vot si6b beviupey srt eras 6p dots 
bevati6p mead 26 stab novos der F089 ofit azth of steb ond esttow 7 7 

and prian 290n6!sd yprens. bis 225m emotiag Moti JWIAS zevsup SATAG 7 
29s 19vs Sldstrev bas 2o0neisd ant ‘etiqcdue has of FF Aath ofid nt si6b oy 
yotoivg snarl oat of nf Jugduo bis s105-¥ oF beveup ets emstpoTg oT 7 
tsded ant vd bennotseq 27 svode bentfiiw nottonu? boesse sit ; 

sonsisa ypyers bis 226m Svidsiss( ne zesob 31 .VADE2 iMsyporg 2299019 . 
bas wort woosy tuqtuo ent .wolt biuprl fugnt sav novrp opsd2 1286 AO al 

mus ni bevioe ad of opst2z nose ewolls 2riT .metstzogmos bruptf ons » & 


noftdt2oq s2fat 5 2eo2u J4D3A entivordue sit .omufos sad To god ond ort 


.efottsups spst2 ons ayer oF tovs92 a 

2in61pos4d Forsnod brs norttplumasA sfsC aoe 

ad [Ttw sobs bine yiopstsea 2tdz nf emarpoiny Eris S16 sont 
ad bfyow yard A>trw nf ysbto Sri} vistemtxdidgs nt . awd at penotina 1 
2queup ylaviditsqey noidw tomit ovpvbred SAT .quy 3299 6 per ben 
yd baqqote brs bedyste 2F ,02MO8..msyporq Jortno> bas ‘notdst y208 ort 
2etds Yd Joe ozI6 ov6 61sb 9F632 Waste bas Toadies. to oad oo | «S2HOA 


‘yd £0MO9 STi} detb and osnt bast ov eso 7sism Forno amt at 

fayosint brs brewotbost ise maior a aut | T2Many 
hee | es We 
5d6b si6t2 sz to sai ldque sis seotetaem fo 
bbs oT _amrtbe a 

7 : 7 7 7 a | 
: a 

Fy - a 


173 


and control matrices stored in RGMO3 can be listed on the printer by 
using the program RGM54. 

As mentioned previously, the accumulation and control action 
is carried out by the program RGM50. The possible options are: 

1) data accumulation only 

2) porportional feedback control 

3) proportional feedback & feedforward control 

4) proportional & integral feedback control 

5) proportional & integral feedback & feedforward control. 
The data is accumulated in the file RGMO3. The data in this file is 
sorted rewritten into disk file RGMO4 by the program RGM46. The data 
in RGMO4 can then be punched onto cards by GM46A, along with a title 
read from RGMO1 and the number of data points stored for each variable. 
The data can then be plotted on the digital plotter for visual inspec- 


tion by using the program RGM53. 


Lb-4 Model Derivation and Reduction Programs 

The terms in the coefficient matrices of the model as pre- 
sented in Chapter III of this work consist of steady state operating 
data. These matrices were calculated from the steady state data py the 
programs DCMOD, DCMD1, and DCMD2. The input required is outlined in 
the program comment statements, and the discrete form of the mode] 
equation is punched onto cards. 

The twenty state model was reduced to a fourth and to a second 
order system. The GEMSCOPE system [54] was used to obtain the fourth 
order model. The second order system was obtained by using the program 
DCRED in which the subroutine LSRED used the least squares method of 


Anderson [ 2 ] to reduce the model. The input data is the original 
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model, the disturbance to be used and the states to be retained. The 
subroutine LSRED was written in general form to allow its adaption to 


any order model and set of disturbances. 


D-5 Miscellaneous Utility Programs 

There are several programs associated with the distillation 
column which don't fall into the above categories. The following is 
a brief description of these programs. 
CRASH/CRAS1: The CRASH system (Computerized Remote Automatic Shutdown) 
is a system of software and hardware which is used to shut down pilot 
plant equipment in emergency situations. It is connected to the dis- 
tillation column and is used as a monitor over periods of operator 
absence. For a complete description of the system, variable limits 
and required input see the CRASH/CRAS1 Dacs Centre System Documentation 
(R.G. McGinnis Nov. 1971). 
RGM43: This program is used in conjunction with the CRASH system and is 
used to set the shutdown ECO on or off. After CRASH has caused a shut- 
down, the ECO is reset by this program, so the pumps can be restarted, 
The program can also be used to simulate the effect of a CRASH shutdown. 
RGM45/CHRTM/RGM47: This system of programs is used to control the opera- 
tion of the bottoms gas chromatograph on the distillation column. It 
initiates sampling, collects the data, and calculates the composition. 
For a listing and description of the programs and hardware and their use 


see the 'Gas Chromatograph Operating Manual' (Revised Nov. 1971). 
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